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Summary — The spectacular successes of the American chemical industry in recent years should not lead us to 
think that the European chemical industry is in any way inferior. On the contrary, working under its 
present conditions, with limited raw materials and energy resources, European industry still maintains a 
good position. In many fields, by the quality of its products or unit cost price, it is the competitor of Ameri- 
can industry. In the sphere of Chemical Engineering, creative inspiration remains always the real source of 
progress. 

If chemical engineering is not yet recognised everywhere in Europe as an independent technology, the 
great number of technologists who contribute every day to the success of chemical industry, are in fact be- 
having as Chemical Engineers. Like Motiéres Bourgeois Gentilhomme, who wrote prose without knowing 
it, they practise a new technology while ignoring its real name. It is always necessary to know something 
of the grammar to speak a language well. The coming generation of technologists enlisted by the industry 
every year, can no longer ignore the principles of a science that will teach them how to conceive, calculate, 
draw, build and operate the equipment by which any kind of chemical reaction may be carried out on an 
industrial scale. 

Wherever it may be—among constructors, in the laboratories, or in factories—original research work 
in the chemical engineering field is constantly being carried out in Europe. **‘Chemical Engineering Science— 
Génie Chimique” is being published to give our colleagues a medium for the publication of their results in 
their own language. 

We hope to publish original work, not only on basic operations of chemical engineering, but also on the 
theoretical principles which rule them: hydrodynamics, thermodynamics ete. We shall also be pleased to 
publish notes on research on industrial processes: constributions which must not be limited to a more or 
less detailed description of the plant, but which must be accompanied by complete manufacturing balances. 

It is to be hoped that the leaders of European chemical industry will give their collaborators every op- 
portunity to publish their work. The problem of “manufacturing secrets” has been the subject of much 
controversy in past years. With the great majority of my university colleagues I, personally, am convinced 
that, in the long run, it is “ot possible to hide one’s light under a bushel; that it is in the end a bad policy 
to keep jealously to oneselt the few grains of truth of which one has temporary possession. 


Le GENIE CHIMIQUE peut étre défini comme 
la branche spéciale de la chimie qui a pour objet de 
foncevoir, calculer, dessiner, faire construire et faire 
conctionner \'appareillage dans lequel on réalisera une 
réaction chimique queleconque a I’échelle industrielle. 
Dans l'art de I’ Ingénieur, le GENTE CHIMIQUE* doit 
venir s'ajouter aux spécialités bien connues du Génie 
Militaire, du Génie Maritime, du Génie Civil ou méme 
du Génie Rural. 

Nous revendiquons pour le GENIE CHIMIQUE, 
le caractére d'une véritable science. Faisant la syn- 
thése de branches particuliéres de la mécanique, de la 
physique ou de la chimie, il apporte au technicien 
l'ensemble des connaissances qui lui sont indispensables 
pour mener a bien la fabrication de immense variété 
des produits industriels dont la préparation est rendue 


* Nous avons été amenés a introduire ce terme nouveau 
de «GENIE CHIMIQUE» pour rendre d'une maniére correcte 
le terme de CCHEMICAL ENGINEERING» devenu, depuis 
bientét trente ans, d'un usage courant dans les pays anglo- 
saxons et spécialement |’ Amérique du Nord. 
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possible par un ensemble plus ou moins complexe de 
réactions chimiques. 

Formant une véritable discipline scientifique avec 
ses méthodes propres pour résoudre ses problémes par- 
ticuliers, le GENIE CHIMIQUE nous vient, en fait 
des Etats-Unis of il a pris naissance et ott son dé- 
veloppement trés poussé explique en partie, sinon 
totalement, les succés remarquables de industrie 
chimique américaine pendant ces derniéres années et 
spécialement durant la guerre. 

Dans la technique européenne, chaque fois qu'il 
s'agit de prendre une réaction mise au point au labo- 
ratoire et de la transposer dans la pratique industrielle, 
on sait qu’il faut de nombreuses étapes et beaucoup de 
travail. On passe d’abord des quelques grammes du 
laboratoire 4 quelques kilogs par jour, puis 4 quelques 
centaines de kilogs, et enfin 4 quelques tonnes. D’une 
maniére générale, dans ce développement progressif 
d'une réaction de laboratoire vers la fabrication in- 
dustrielle, on procéde par sauts successifs, avee un 
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facteur de transposition voisin de 10: au grand maxi- 
mum, on passe du kilog 4 heure A la tonne par jour, 
soit un facteur de 40. Quand on a un atelier qui fait 
une tonne par jour, on se lance alors a faire A un 
atelier dix fois plus grand, et l'on juxtapose deux 
ateliers, trois ateliers, pour faire face aux besoins de 
la production. 

Les dures nécessités de la deuxiéme guerre mondiale 
ont obligé lindustrie chimique américaine a éxécuter, 
en quelques mois seulement, cette transposition d'une 
réaction de laboratoire vers la fabrication industrielle. 
Production colossale de caoutchouc synthétique a par- 
tir du pétrole, fabrication intensive de produits entiére- 
ment nouveaux (D.D.T., pénicilline, ete.), tout cela 
a di étre monté de toutes piéces en moins de deux ans. 
Que dire de la réalisation de la bombe atomique ex- 
igeant la fabrication industrielle d’isotopes purs: la 
premiére pile atomique dégageant 1200 kWh était en 
marche moins d'un an aprés que la premiére pile de 
laboratoire, mettant en jeu une puissance de 0.2 kWh, 
ait produit les quelques microgrammes de plutonium 
indispensable aux caleuls. Pour obtenir de tels ré- 
sultats, il a fallu que l Amérique disposat de méthodes 
de travail, d’états-majors, de cadres de techniciens 
capables de concevoir, de calculer, de dessiner, de faire 
fonctionner toutes ces installations nouvelles. 

Si durant la deuxiéme guerre mondiale, |’industrie 
chimique américaine a pu surclasser nettement l’in- 
dustrie chimique allemande, nous affirmons que c'est 
grace A l’avance prise par les Etats-Unis dans le do- 
maine du GENIE CHIMIQUE. Pourtant nos ad- 
versaires disposaient certainement d’états-majors et 
de cadres de haute qualité, mais dans cette bataille 
d’états-majors et de cadres, c’est l'emploi de la bonne 
méthode de travail, de la doctrine technique éprouvée 
qui a décidé de la victoire des Américains. 

Aux Etats-Unis, le GENIE CHIMIQUE existe 
comme discipline spécialisée depuis une trentaine d'an- 
nées. Sa doctrine et ses méthodes scientifiques sont 
enseignées dans toutes les Universités importantes. Le 
GENIE CHIMIQUE n’a pas regu ses lettres de na- 
turalisation en Allemagne en tant que discipline scienti- 
fique séparée. Quant A la France ot |'Université de 
Toulouse a été la premiére & consacrer le terme du 
GENIE CHIMIQUE, la signification de ce titre n'est 
pas toujours saisie par les industriels et par les chi- 
mistes. 

Il n’est pas mauvais d’examiner le concours des 
circonstances qui ont amené les Etats-Unis 4 prendre 
en si peu de temps, une avance si considérable et un 
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retour sur le développement de l'industrie chimique 
moderne permettra de saisir la période cruciale oi dans 
les deux continents, les techniques industrielles du chi- 
miste ont pris des routes divergentes. 

Au XIX éme siécle, tandis que l'industrie mé- 
canique évoluait trés rapidement vers une technique 
raffinée grace A la thermodynamique et aux développe- 
ments tirés du principe de Carnot, l'industrie chimique 
est restée un art, bien qu’d la suite de Lavoisier, Vau- 
quelin, Berthollet, ete., elle ait quitté le domaine de 
Vartisanat dés le début du 19 éme siécle. C'est simple- 
ment que devant créer une industrie 4 partir de ré- 
actions étudiées au laboratoire, le chimiste ne disposait 
pas de matériaux de construction ou d'outils de tra- 
vail, basiquement différents de ceux qu'il utilisait au 
laboratoire. La céramique ou la magonnerie venaient 
remplacer le verre, le tirage d'une cheminée faisait 
circuler les gaz: les liquides corrosifs circulaient par 
gravité, plus tard par monte-jus, lorsque les progrés de 
la mécanique mirent l'air comprimé a la disposition du 
chimiste. 

Par un effort tenace, avec énormément d’ingénio- 
sité, les chimistes arrivérent ainsi 4 mettre sur pied des 
industries presque parfaites: procédé des chambres de 
plomb, procédé Solvay, etc. qui subsistent encore prati- 
quement inchangés aprés un siécle ou presque. Dans 
ces procédés, le tour de main, le secret de fabrication 
restent encore la marque typique du parfait praticien 
rompu par une longue pratique a tous les aléas du 
métier. 

C'est lavénement du «procédé de contact» pour la 
fabrication de acide sulfurique qui marque, croyons- 
nous l’avénement de l'industrie chimique moderne. La 
fabrication de l’indigo synthétique exigeait la produc- 
tion a bas prix de loléum, obtenu jusqu’alors par le 
vieux procédé de Nordhausen. Les travaux de 
Knretscu, pour le compte de la «Badische Anilin und 
Soda Fabrik» aboutissent en 1900 a la fabrication in- 
dustrielle de 'anhydride sulfurique par oxydation cata- 
lytique directe du gaz sulfureux, suivant une réaction 
étudiée par KUHLMANN vers 1830. 

Au temps de KuBLMANN, l'art chimique n’était pas 
assez évolué pour permettre de réaliser la purification 
a lextréme des gaz de grillage des pyrites. Il était 
également impossible de réaliser un contréle précis 
de la température des appareils de synthése. _KN1ETSCH 
pouvait disposer d’appareils précis, sensibles, écono- 
miques pour mesurer la température de ses gaz: le 
couple thermoélectrique de Le CHATELIER avait fait 
ses preuves. La technique allemande, déja trés en 
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avance dans la synthése organique, s'est attaquée et a 
vaincu les difficultés de la purification des gaz. Jusqu’a 
la premiére guerre mondiale, elle garda le monopole 
presqu’exclusif de la fabrication de loléum pour la- 
quelle elle monta des usines tres considérables, pro- 
duisant des dizaines de tonnes par jour. 

Essentiel pour la réussite de la fabrication, le con- 
trole rigoureux de la température exige immédiatement 
que les gaz circulent avec des vitesses bien constantes, 
aisément réglables dans toute linstallation. Il faut 
également faire circuler dans des conditions parfaite- 
ment contrélables, l'acide destiné a fixer l'anhydride 
sulfurique formé. Vannes, pompes, ventilateurs de- 
viennent les appareils essentiels dont dépend tout le 
succés. Le chimiste doit faire appel au mécanicien qui 
devient un personnage important, car tout dépend 
de la bonne marche des installations mécaniques, de 
leur entretien malgré lusure ou la corrosion. 

Durant la premiére guerre mondiale, l'industrie 
chimique allemande fut appelée A développet encore 
sa collaboration avec les ingénieurs mécaniciens. Pour 
industrialiser le procédé HABER pour la fabrication de 
lammoniaque synthétique, il fallut construire des 
compresseurs géants capables de fonctionner sans arrét 
pendant des centaines d’heures pour comprimer a 
250 atmosphéres les milliers de tonnes de gaz 
4 transformer. L’industrie chimique dut faire cons- 
truire de véritables canons, mais des canons dont les 
dimensions dépassaient de beaucoup celles des 420 ou 
de la «grosse Bertha». A la demande du chimiste, le 
mécanicien et le métallurgiste résolurent trés rapide- 
ment tous ces problémes. 

En Allemagne, de 19144 1918, la collaboration entre 
mécanicien et chimiste avait conduit a l'ammoniaque 
synthétique, A Vacide nitrique synthétique faisant 
appel a des alliages nouveaux résistant a la corrosion: 
les aciers inoxydables. Aprés sa défaite, Allemagne 
s'est laneée dans la fabrication des carburants synthé- 
tiques remplagant le pétrole par hydrogénation de la 
houille, ensuite dans la fabrication du caoutchouc syn- 
thétique, puis dans la fabrication des matiéres grasses 
encore A partir du charbon par le procédé OXO, ete. 

Durant la période 1920-1944, le travail du chi- 
miste allemand a été un travail absolument remar- 
quable, car il a abouti A créer toute une série nouvelle 
de réactions de laboratoire. La collaboration entre 
chimiste et mécanicien est devenue de plus en plus 
féconde et a permis de transposer toutes ces réactions 
dans la pratique industrielle, de sorte qu’A l'heure actu- 
elle, on peut dire sans exagération que l'industrie chi- 
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mique peut faire n’‘importe quel produit avec n’im- 
porte quelle matiére premiére, si elle dispose d’éner 
gie et de main-d’oeuvre en quantité suffisante, lorsque 
la question du prix de revient n’intervient pas. 

Durant la guerre de 1914—1918, les chimistes fran- 
cais et britanniques réussirent par un effort admirable 
a créer de toutes piéces lindustrie chimique indispen- 
sable & la victoire. Ici aussi, il fallut faire un appel de 
plus en plus fréquent au concours d’ingénieurs mé- 
caniciens, Ces derniers, malheureusement, ont trés rapi- 
dement pris la prépondérance sur le chimiste, aussi 
bien en France qu’en Grande-Bretagne, tandis qu’en 
Allemagne, chimistes et mécaniciens collaboraient sur 
un pied d’égalité. Mais, en Allemagne, comme en France 
et comme en Grande-Bretagne, il y a deux branches 
distinctes de techniciens, de formation trés différente. 
Trés souvent, elles ont collaboré trés efficacement, 
c'est le cas de l’Allemagne, quelquefois, l'une a écrasé 
Vautre, c'est le cas de la France et de la Grande- 
Bretagne . 

C’est A cause de cette dualité que la discipline dont 
nous nous oceupons ici, le «GENTE CHIMIQUE), 
ne s'est pas créée en EUROPE. C’est par un concours 
plus heureux des circonstances que cette discipline 
s'est créée aux Etats-Unis. 

Sur le plan mondial, en 1914, la chimie scientifique 
ou industrielle des Etats-Unis ne comptait que pour 
fort peu. Dans les techniques nouvelles de l’automobile 
ou de aviation, comme pour la création de leur in- 
dustrie chimique de guerre, nos alliés durent faire 
appel A nos techniciens, copier les procédés montés 
en France ou en Grande-Bretagne pour créer leurs 
industries d’armement. 

De 1920 A 1945, la situation s'est complétement 
renversée et sitét la seconde défaite de Allemagne, 
c'est en sens inverse que sont allées les missions d’étu- 
des, en quéte de techniques nouvelles, achetant fort 
cher les licences de fabrication, important des «know- 
how», en méme temps que le «slang» américain. 

Pendant la guerre 1914—1918, il a fallu, aux 
Etats-Unis comme en Europe, mettre en commun 
toutes les ressources, faire appel A toutes les bonnes 
volontés, faire collaborer professeurs, mécaniciens, 
chimistes de laboratoire, chimistes industriels. Mais 
sous influence de quelques personnalités trés mar- 
quées, cette collaboration entre chimistes et mécani- 
ciens s'est établie sur des bases toutes différentes. 
elle a trés vite abouti 4 une fusion. Le «CHEMI- 
CAL ENGINEERING», le «GENIE CHIMIQUE» 
est devenu une discipline scientifique reconnue. La 
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profession de «CHEMICAL ENGINEER», dingénieur 
du GENTE CHIMIQUE a acquis son domaine distinct 
dans les spécialités diverses du métier d’ingénieur. 
C'est par une remarque trés astucieuse que nos coll- 
égues américains ont pu faire d’énormes progrés dans 
’établissement sur des bases scientifiques de la techno- 
logie chimique lorsqu’ils ont remarqué que toute fabri- 
cation chimique, quelles que soient les matiéres pre- 
miéres qu'elle traite, fait passer ces diverses matiéres 
par une succession d’opérations distinctes. Ces opé- 
rations restent toujours identiques a elles-mémes dans 


leur principe si les modalités d’applications varient 
suivant la nature des produits fabriqués. C'est de la 
qu'ils ont tiré le concept véritablement original de 


«Unit operation» que nous traduirons, faute de 
mieux, par «Opération Fondamentale». 

Il semble aujourd’hui évident de penser que la sé- 
paration d'un liquide d’aveec les matiéres solides qu'il 
contient en suspension (filtration) constitue une opé- 
ration proprement dite, qui doit étre étudiée en soi, 
tout A fait indépendamment de l'industrie chimique 
particuliére dans laquelle elle s applique, qu'il sagisse 
par exemple de la clarification des eaux potables, de 
la séparation des cristaux au sein d’une eau mere, ou 
de la filtration d'une pulpe de flottation. Mais il a 
fallu pourtant plus d'un siécle pour que les techniciens 
de Vindustrie chimique échappent A lempirisme de 
leur propre spécialité, définie par le produit chimique 
quils devaient personnellement fabriquer (acide sul- 
furique, sucre, extrait tannant ou métal) pour établir 
les lois fondamentales qui réglent un tel phénomeéne. 
A notre avis I'«Unit Operation» ce concept si simple 
que sa découverte peut paraitre maintenant un véri- 
table «Oeuf de Colomb», marque pour l'industrie 
chimique la rupture compléte d’avec les procédés de 
«homme de l'art», il fait rentrer toutes les opérations 
si complexes des fabrications chimiques dans le cadre 
solide de la technique scientifique. 

Par «Unit Operation», il faut entendre une opé- 
ration industrielle essentiellement indépendante. Cette 
opération obéira A des lois physiques ou chimiques par- 
faitement définies et se retrouvant inchangées dans 
un grand nombre de procédés de fabrication divers. 
C'est surtout par quelques exemples que l'on pourra 
mieux comprendre ce que les Américains entendent par 
«Unit Operation» puisque la langue frangaise ne pos- 
séde pas, A notre connaissance, de terme qui rende 
parfaitement bien ce que nos alliés veulent entendre 
par la. «Operation» pour nos alliés, cela indique d'une 
fagon trés générale, une série de manipulations, mais 
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le terme «Unit» vient y ajouter Vidée que cette série 
de manipulations forment un tout, un ensemble qui 
posséde en propre une identité particuliére. Le terme 
«Operation fondamentale» nous parait rendre a peu 
prés le méme sens. L’industrie chimique doit reposer 
sur l'emploi systématique d’opérations fondamentales 
trés diverses, tant physiques que chimiques. 

Par exemple, dans un trés grand nombre d’opé- 
rations industrielles, les mati¢res premiéres devront 
étre divisées pour augmenter leur surface: cette opé- 
ration, dans le cas des solides, c'est le broyage qui 
forme une classe particuliére d’«Unit Operation» pour 
le GENTE CHIMIQUE, S’il s’agit de matiéres liqui- 
des ou gazeuses, cette augmentation de la surface 
sobtiendra soit par émulsification, soit par atomisa- 
tion, autres types d’opérations fondamentales. II est 
de la plus haute importance pour l'Ingénieur du GE- 
NIE CHIMIQUE, responsable d'une fabrication in- 
dustrielle donnée, de connaitre les lois théoriques et 
pratiques qui régissent les diverses ¢tapes intermé- 
diaires du procédé dont il a la responsabilité. 

Parmi les «Unit Operations» les plus importants. 
nous indiquerons (’abord les plus simples qui se con- 
tentent de séparer diverses substances sans change- 
ment de phase: séparation électrique ou magnétique. 
séparation hydraulique, flottation, décantation, fil- 
tration, centrifugation ete. Des «Unit Operations» 
plus complexes comporteront des séparations entrai- 
nant en méme temps, un changement de phase pour les 
diverses matiéres premiéres: dissolution ou extraction, 
évaporation, distillation, séchage, cristallisation, ab- 
sorption, conditionnement de lair, ete. 

D’autres Opérations Fondamentales sont caracté- 
risées par Videntité du processus purement chimique 
auquel il est fait appel: hydrogénation, oxydation, 
nitration, sulfonation, estérification des composés de 
carbone par exemple *. 

Mais pour le GENIE CHIMIQUE, la distinetion 
entre Chimie Minérale ou Chimie Organique ne con- 
serve aucun sens. La synthése de 'ammoniac et celle 
du méthanol font appel aux mémes opérations fonda- 
mentales et s'effectuent dans des appareillages indus- 
triels de méme type. L’oxydation du méthanol en 
formaldehyde, celle de lammoniac en oxyde nitrique 

* Les américains semblent maintenant réserver le nom 


@eUNIT Process» aux opérations fondammentales mettant en 


jeu des transformations chimiques. lls conservent le nom 
d’«Unit Operation» aux opérations fondamentales compor- 
tant exclusivement des transformations physiques. A lorigine, 
la distinction entre «Unit Operations et «Unit Processs était 


loin d'étre universellement adoptée. 
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appartiennent au méme type d’Operation Fonda- 
mentale. 


Au point de vue pédagogique, il est certainement 
beaucoup plus logique et beaucoup plus fructueux 
d’étudier ’Opération Fondamentale ou «Unit Opéra- 
tion» quest la distillation comme un procédé géné- 
ral de transformation, plutét que de létudier dans le 
cas particulier de la rectification des jus alcooliques ou 
de la séparation du goudron qui intéressent plus parti- 
culiérement les ingénieurs des industries de fermenta- 
tion ou ceux de Vindustrie gaziére. Si certaines Opéra- 
tions Fondamentales, comme [agitation ou le brovage 
sortent A peine de lVempirisme, d'autres au con- 
traire, reposent sur des bases théoriques parfaitement 
connues aujourd’hui. Pour les Opérations Fondamen- 
tales de l’évaporation, de la vaporisation, de la dis- 
tillation, de absorption ou de extraction, le GENIE 
CHIMIQUE fournit A lVingénieur des méthodes de 
caleuls dont la précision est parfaitement connue et 
parce qu'elle dispose d’Ingénieurs nombreux pos- 
sédant la maitrise compléte des méthodes du GENIE 
CHIMIQUE que lindustrie chimique américaine a pu 
surclasser indiscutablement lindustrie chimique alle- 
mande. 

Mais, ce serait méconnaitre gravement le caractére 
du GENIE CHIMIQUE si on ne le considérait que 
comme une branche un peu spéciale (et probablement 
inférieure) de la Physico-Chimie, sous prétexte qu’il 
connait ses succés les plus nets dans des opérations 
fondamentales dont les lois thermodynamiques sont 
parfaitement connues. 


Bien au contraire, en étudiant individuellement 
chacune de ses opérations fondamentales sous langle 
qui lui est particulier, et avec les méthodes qui lui sont 
propres, le GENTE CHIMIQUE arrive A des géné- 
ralisations fécondes. Dans toutes ces opérations de 
distillation, extraction, dissolution ete., il trouve un 
exemple particulier des «ephénoménes d’échange» dont 
l'étude et la théorie me paraissent lui appartenir en 
propre. 

Dans ces phénoménes dont la cinétique surtout lui 
importe, il envisage le passage d’espéces moléculaires 
différentes au travers d'une surface de séparation. En 
cherchant a utiliser les «coefficients de diffusion» sur 
lesquels la physique théorique ou expérimentale lui 
fournissent si peu d’informations réelles, il est con- 
duit a utiliser en les généralisant d'autres «ephénoménes 
d’échange» appartenant a des domaines propres de la 
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dant, ne peut étre réduit 4 un simple compartiment de 
la Physique Appliquée. 

Ce nouveau concept de «phénoménes d’échanges» 
s'introduit tout naturellement lorsqu’on cesse d’ex- 
aminer l'industrie chimique en général du point de vue 
analytique qui nous a permis de décomposer chaque 
fabrication en «Opérations Fondamentales». Nous 
pouvons également considérer l'industrie chimique d’un 
point de vue synthétique et chercher ce qui est abso- 
lument général dans toute espéce d’industrie chimique. 
L’on peut affirmer sans crainte de contredit que toute 
fabrication chimique est d’abord un probléme de circu- 
lation: toute fabrication doit résoudre des problémes 
de transport. 

Dans toute fabrication industrielle, il s’agira iné- 
vitablement d’introduire, puis de retirer d'un appa- 
reillage approprié les matiéres rentrant dans le cycle 
de transformation. Dans le cas particulier des fluides, 
liquides ou gaz, le GENIE CHIMIQUE est done ob- 
ligé de connaitre 4 fond les méthodes scientifiques de la 
Mécanique des Fluides. En particulier, il est obligé 
d'utiliser les relations de simi):tude et les méthodes 
d’approche de l’analyse dimentionnelle avec l'utili- 
sation de coefficients sans dimensions, tel que le nom- 
bre d’Osborne Reynolds. 

Dans toute réaction chimique, il est nécessaire égale- 
ment, soit d’apporter, soit de soustraire a l'appareillage 
réactionnel une certaine quantité de chaleur. Les 
phénomeénes d’échange thermique, soit par rayonne- 
ment, mais surtout par convection, appartiennent 
également au domaine du GENIE CHIMIQUE. 

Dans ces deux genres de problémes d’échange, on 
retrouve des équations de forme analytique semblable : 
les relations de similitude introduisent ici les nombres 
sans dimension de de ou de Gras- 
Hor. Pour rendre compte des phénoménes de convec- 
tion, on est conduit A introduire le concept d’une 
barriére particuliére, d'une couche de passage, d’un 
film créant une discontinuité a la surface de séparation, 
film que les expériences sur |’écoulement des fluides 
permettent de rendre visible. 

Pour interpréter la diffusion d’espéces moléculaires 
diverses au travers de la surface de séparation existant 
entre deux phases, le GENIE CHIMIQUE, est tout 
naturellement amené a utiliser encore les propriétés de 
ce film particulier, barriére naturelle qui, pour lui 
s’oppose aussi bien aux échanges de quantité de mou- 
vement qu’aux échanges de qualité de matiére: a 
l'échelle de la molécule, en fait, tous ces phénoménes 
d’échange dépendent des contacts individuels entre 
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particules réelles. Le frottement éprouvé par un fluide 
s’écoulant au contact d'une paroi solide résulte en fait 
des perturbations subies individuellement par chaque 
molécule heurtant la paroi. Au cours de choes plus 
ou moins parfaitement élastiques, sa quantité de mou- 
vement varie et le résultat statistique de ces variations 
perturbe la composante moyenne de son déplacement 
par rapport A la paroi. Dans léchange thermique, le 
processus ¢lémentaire reste le méme: dans le choc, 
lénergie cinétique d'une molécule augmente ou dimi- 
nue au contact des molécules plus lentes ou plus ra- 
pides, plus chaudes ou plus froides du métal constitu- 
ant la paroi. Dans les échanges entre vapeurs et liqui- 
des complexes (si importants pour la distillation et la 
rectification), les passages individuels plus nombreux 
d’une espéce moléculaire au travers de la surface de 
séparation entre les deux phases déterminent statisti- 
quement dans chaque phase les variations de concen- 
tration qui intéressent le chimiste. 

Par une méthode scientifique qui lui est propre, le 
GENTE CHIMIQUE arrive 4 interpréter de la méme 
fagon, par les mémes opérations mathématiques, les 
phénoménes en apparence trés divers que le méme con- 
cept physique lui a permis de rapprocher. I! prend 
done rang de véritable discipline scientifique, faisant 
sa synthése particuliére de chapitres tout A fait di- 
stincts de la physique ou de le physico-chimie. 

Basant sa formation professionnelle sur une solide 
culture mathématique, physique et chimique, l'in- 
génieur du GENTE CHIMIQUE doit acquérir une men- 
talité particuliére, posséder des réflexes intellectuels 
spéciaux qui lui permettront d’aborder avec fruit les 
problémes caractéristiques de sa profession: le facteur 
temps constitue un des éléments essentiels de ses 
problémes. 


A la différence du thermodynamicien qui se pré- 
occupe avant tout de léquilibre qui viendra figer toute 
transformation, lingénieur du GENIE CHIMIQUE 
devra reprendre 4 son compte les préoccupations du 
cinéticien qui va chercher a déterminer le temps né- 
cessaire pour lachévement de cet équilibre. Signalons 
immédiatement que la cinétique propre au GENIE 
CHIMIQUE est beaucoup plus complexe que la ciné- 
tique classique. 

Tandis que le pur cinéticien peut détailler 4 volonté 
le mécanisme réactionnel en en fixant étroitement tou- 
tes les conditions (température, pression, concentra- 
tion, ete.) lingénieur du GENTECHIMIQUE se trouve 
devant une tache beaucoup plus complexe: il ne 
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peut rester maitre des conditions opératoires, il ne 
peut plus traiter individuellement chacun des phéno- 
ménes d’échange. Dans le probléme réel, il est impos- 
sible de les dissocier les uns des autres. L’élimination 
des calories produite par la réaction commande |’in- 
fluence que la température prendra sur la vitesse de 
réaction; cette température agit aussi sur les coeffi- 
cients de diffusion et réagira obligatoirement sur la 
concentration des molécules actives A la surface du 
catalyseur. Le probléme concret du caleul et du 
dessin convertisseur catalytique ne peut ¢étre 
mené da bien qu’en tenant un compte rigoureux de 
toutes ces influences diverses, réagissant les unes sur 
les autres. 

Parfaitement outillé maintenant par la connais- 
sance de tous les détails qui conditionnent chaque Ope- 
ration Fondamentale, lingénieur du GENTE CHIMI- 
QUE doit encore intégrer l'ensemble de ces données 
pour établir le «evele de fabrication» qui dirigera dans 
un ordre rigoureux chacun de ces plans de détail. Il 
saura d’avance le retentissement de chaque variable 
sur le résultat d'une opération de détail, il connaitra 
celles qu'il faut controler rigoureusement et se sera 
préoceupé des moyens indispensables pour assurer ce 
controle. 

La création d'un cycle de fabrication est en fait un 
véritable travail d’architecte. Dans une Ecole des 
Beaux-Arts, on n’apprend pas aux jeunes architectes 
& faire A coup sir des chefs-d’oeuvre, mais on leur 
donne des principes sirs qui leur permettront de faire 
des constructions solides, adaptées 4 leur objet, La 
maitrise de Varchitecte déterminera lélégance de 
ladaptation de la construction A son objet. Le GE- 
NIE CHIMIQUE fournit des régles sires pour ce tra- 
vail essentiel. 

L’architecte travaille sur plans, Vingénieur du 
GENIE CHIMIQUE a des plans particuliers 4 sa 
technique: ce sont les bilans. On ne peut pas conce- 
voir un cycle de fabrication sans faire immédiatement 
le bilan des opérations, bilan qui doit établir ce que 
devient toute matiére introduite dans la suite des 
étapes de transformation. Le bilan mati¢re doit rappeler 
4 Vingénieur du GENIE CHIMIQUE que la loi de 
LAVOIsIER est indiscutable. Une réaction chimique 
est rarement parfaite: la plupart du temps elle est ac- 
compagnée de réactions parasites qui viennent dimi- 
nuer le rendement. Le bilan matiére doit tenir un 
compte exact de tout ce qui est détourné A une étape 
queleonque de la fabrication. En fait Vingénieur du 
GENTE CHIMIQUE doit étre un comptable, mettant 
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constamment en balance ce qui entre et ce qui sort 
d'un atelier. A cété des pertes occasionnées par des 
réactions parasites, il faut tenir un compte exact des 
pertes dues 4 des opérations imparfaites: lavage en- 
trainant une certaine fraction de produits solubles, 
pertes en queue pour des gaz imparfaitement absorbés, 
etc. également aux imperfections de l'appareillage qui 
présente des fuites. 

A cété du bilan matiére, il y a un autre bilan tout 
aussi essentiel, c'est le bilan énergie: Vénergie coite 
cher comme les mati¢res premiéres. Ce bilan énergie 
se fait exactement comme le bilan matiére, en caleu- 
lant A chaque étape les quantités de chaleurs intro- 
duites avec les diverses matiéres premiéres, ou les 
divers moyens de chauffage et en retrouvant sous 
forme de produits chimiques transformés les quantités 
d'énergie introduites. 

Cette notion du prix de chaque chose vient rappeler 
opportunément que Vingénieur du GENIE CHIMI- 
QUE ne peut se contenter d’étre un pur théoricien: 
il doit constamment rester pénétré du souci du pos- 
sible. Le prix de revient doit toujours demeurer sa 
préocecupation dominante, et ce souci ne peut que mo- 
difier profondément sa mentalité propre par rapport a 
celle du chercheur, du physicien ou du chimiste de la- 
boratoire. 

Dans l'industrie chimique, nous ne prétendons pas 
que lVingénieur du GENIE CHIMIQUE doive étre 
partout présent et exclusivement présent, mais nous 
affirmons qu'il doit lui étre fait une place importante a 
cété du chercheur dans son laboratoire de recherches. 
de Vélectricien ou du mécanicien dans la centrale ou 
latelicer d’entretien qui leur sont confiés, A cété de 
l'économiste qui surveille les marchés et la production 
de lensemble industriel, A c6té de ladministrateur qui 
décide la politique de organisation et doit prévoir son 
avenir. L’ingénieur du GENIE CHIMIQUE a droit 
Ad une place dans les bureaux d’étude, dans les ateliers 
de fabrication et dans les états-majors del industrie 
chimique. 


Les réalisations spectaculaires de l'industrie chimi- 


que américaine pendant ces derniéres années ne doi- 


vent pas faire croire que Vindustrie chimique euro- 
péenne est en pleine décadence. Bien au contraire, 
travaillant dans les conditions qui lui sont propres, 
avec les matiéres premiéres et les ressources en énergie 
dont elle dispose, elle fait encore trés bonne figure. 
Dans bien des domaines, soit par la qualité des pro- 
duits fabriqués, soit par le prix de revient unitaire, 
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elle concurrence victorieusement l'industrie d’Outre- 
Atlantique. Dans l'industrie chimique, l'esprit créa- 
teur reste toujours la véritable source du progrés. 


Si le GENTE CHIMIQUE nest pas encore partout 
‘reconnu en Europe comme une discipline scientifique 
indépendante, les trés nombreux techniciens qui con- 
courent chaque jour aux succés de l'industrie chimique 
se comportent en fait en ingénieurs du GENIE CHI- 
MIQUE. Comme Mr. Jourparn faisait de la prose 
sans le savoir, ils pratiquent une discipline scientifique 
nouvelle sans lui accorder son vrai nom. I] est utile de 
connaitre la grammaire pour bien parler une langue. 
Les nouvelles générations de techniciens que lindu- 
strie doit recruter chaque anné ne pourront plus ignorer 
les principes de la science qui leur apprendra A con- 
cevoir, caleuler, dessiner, faire construire et faire fone- 
tionner l'apparaillage dans lequel s‘effectuera a l'échelle 
industrielle, une réaction chimique quelconque. 


Que ce soit chez les constructeurs, dans les labo- 
ratoires industriels, ou dans les ateliers de fabrication. 
il se poursuit constamment, en Europe des travaux 
originaux de recherches dans le domaine du GENIE 
CHIMIQUE. «Chemical Engineering Science —- Ge- 
nie Chimique» va paraitre pour fournir a nos nombreux 
collégues un organe d’expression qui pourra publier 
leurs travaux, dans leur propre langue. 


Nous espérons publier des travaux originaux por- 
tant non seulement sur les Opérations Fondamentales 
du GENTE CHIMIQUE, mais également sur les prin- 
cipes théoriques qui les régissent: mécanique des 
fluides, thermodynamique, etc. Nous serons également 
heureux de publier des études sur des procédés indu- 
striels, études qui ne se borneront pas simplement a 
une description plus ou moins détaillée de l’appa- 
reillage, mais qui seront accompagnées de bilans com- 
plets de fabrications. 


Nous souhaitons que les dirigeants de l'industrie 
chimique européenne donnent a leurs collaborateurs 
toutes les facilités nécessaires pour la publication de 
leurs travaux. La question du «secret des recherches» 
a fait couler beaucoup d’encre ces derniéres années. 
Comme la grande majorité de mes collégues universi- 
taires, je suis personnellement convaincu qu’A la 
longue, on ne réussit jamais A tenir la lumiére sous le 
boisseau et que c’est finalement une mauvaise politique 
que dessayer de conserver jalousement pour soi les 
quelques parcelles de vérité que lon est provisoire- 
ment seul A détenir. 
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Rapid calculation of plate number and reflux ratio in batch distillation 
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Summary—The method of Bowmay-Cicurii: for determination of minimum reflux and plate require- 
ments in binary batch distillations has been extended with a relation between these minimum requirements 
and the number of plates and reflux ratio to be used in practice. In the case that the concentrations in the 
column are in equilibrium with the bottom concentration at any moment during the distillation the effect 
of holdup, too, is accounted for. 

The application of the method of calculation has been simplified by the construction of three alignment 
charts with three additional diagrams. 

An experimental verification showed the method also to be applicable to systems with more than 
two components, if certain conditions are fulfilled. 

The concentrations in the column not in equilibrium with the bottom concentration give rise to an 
actual reflux ratio lower than that caleulated: in some cases it may even be lower than the minimum reflux 


ratio. 


Zusammentassung— Die Bowmas-Cronetot Methode fiir die Bestimmung der Minimum Ricklauf- und 
Bodenanspriiche in’ bindren diskontinuierlichen Distillationsverfahren ist durch eine Bezichung zwischen 
diesen Minimum Anspriichen und der in der Praxis anzuwendenden Bodenzahl und Riicklaufverhaltnis 
erweitert worden. Falls die Konzentrationen in der Kolonne in jedem Augenblick mit der Blasenkonzen- 


tration im Gleichgewicht stehen ist auch der Effekt des ,.Holdups* 


beriicksicht igt. 


Die Anwendung der Berechnungsmethode ist durch die Konstruktion von drei Nomogrammen mit 


drei Zusatzdiagrammen vereinfacht worden. 


Es konnte experimentell gezeigt werden, dal die Methode auch anzuwenden ist fur Systeme mit mehr 


als zwei Komponenten wenn gewissen Bedingungen entsprochen worden ist. 

Wenn die Konzentrationen in der Kolonne nicht mit der Blasenkonzentration im Gleichgewicht stehen, 
hat dies zur Folge, da® das wirkliche Riicklaufverhaltnis niedriger als das berechnete ist; in einigen Fallen 
kann es sogar niedriger sein als dae Minimum Ricklaufverhaltnis. 


INTRODUCTION 


Calculation methods for binary separations by batch 
rectification have been given by various authors. 
Rose and co-workers [9, 10, 12, 13, 14] have shown 
how the distillation curve may be computed if reflux 
ratio, plate number and holdup are known. Bocarr|2}, 
[6, 7) and Yu Cars Cue [17] 
gave methods for calculating the yield of the volatile 
component when the reflux ratio is continuously 
increased in order to keep the purity of the product 
at a specified value. 

Recently, Bowman and developed 
an ingenious method by which the minimum reflux 
and plate requirements in batch distillation may be 
found very quickly. They defined the difficulty of 
separation by a separating factor, called “pole 
height,’ which is obtained by multiplying the slope 
of the distillation curve at the point where the 
distillate contains 50 mol.-°e volatile component by 
the mol.-° of the residue left in the reboiler at the 
same moment: . dy 


=W aw 


(1) 


The pole height « is related by the following simple 
equations to minimum plate number and minimum 


reflux ratio: 
log S« 


N 
log x 


(2a) 


R (2b) 


» 
abt 


The equations are only applic. 
able when holdup of the column is negligible; further 
the relation between the minimum requirements and 
plate number and reflux ratio and the actual values 
to be used is not yet known. 

The present investigation was initiated to develop 
this relation, taking the influence of holdup into 
account*. 


DEVELOPMENT OF THE CORRELATION BETWEEN R, 
R,,, N N,, 

Calculation of batch distillation processes when holdup 

is appreciable is complicated by the fact that the 


* At the same time as this paper was submitted for 
publication, M. T. and co-workers (Ind. Eng. 
Chem. 1950 42 4502-10) published a similar study, in which 
the relation between pole height, reflux ratio and plate 
number is given. Their results, though given in a different 
form, are in agreement with the present work. However, 
Crene tii et al. did not consider the influence of holdup, 
nor Was any experimental work included. 
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concentrations in the column cannot always be con- eq. (5a), (5c) and (6b) are reduced to: 


sidered in equilibrium. This is caused by the con- 


tinuous removal of the most volatile components at r=} - E° (7a) 
the top, resulting in a continuous decrease in concen- 

tration of these components in reboiler and column. dy (1+? (7b) 
Only if certain conditions are fulfilled can the concen- iad " 

trations in the column be interrelated by the well- dm E?—1 

known operating line as used in continuous distillation. dx 2EinE’ (8) 


This fulfilment is obtained when: 
1. The rate of change in concentration 
is small (as in the case of a bad separa- 


Table 1. Calculated pole height values as a function of holdup quotient 
and log E; column at total reflux 


tion or with low rate of product removal), — jog HW 
2. the rate of equilibration of the A 0 0-05 010 O15 0.20 0.25 0-30 0-40 0-50 
column is great (at high liquid and va- ; 
pour loads, or low plate number). O70 0-600 0-522 0-463 0-415 0-376 0-343 
Under these conditions the following 100) -1-238 1-016 0-868 0-752 0665 0-597 
material balance, set up for the case in 1300 2-50 1-87 1-50 - 1250 1072 0-936 


235 206 185 1-67 1-395 1-200 


which a differential amount, of distillate 15 395 207) 271 
200 12-50 760 5-98 475 3-04 3°35, 2-93 2-34 1-94 


is removed at the top of the column, 


16-7 1955 | 7-70 6-09 - 4.27 3-29 2-68 
1055 805 655 5-50 4:17 3-46 

28-8 192 144 9-57 7-20 575 

28-8 23-0 — | 192 14-4 11-50 

- | 288 21-6 17-25 


may be used as a starting point: 2-50 | 30-50 
3000 — 
d(Wr) + Hdx, = ydW (3a) 5-00 
from which wernt 
15-00 
(3b) 
Introduction of the separating factors — W + in (3b) 
results in 
(y 2) dy 
dx (4) 
H dx, 
1+... 


(1) Minimum plate number 

Under conditions of minimum plate number (e.g. 
total reflux), the following relations between top 
composition and bottom composition, and holdup 
composition and bottom coniposition may be used: 


(5a) 
in which 
E=2* (5b) 
and [14] 
Infl + (2 
InE 


From (5a) and (6a) follows 
dy [y+(l—B}* 


E ( 


and 
dx, E—1 


dz (1+ (6b) 


Using the stipulation that the criterion of separa- 
tion is applied when distillate composition y= 0-5, 
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Substitution of these values in eq. (4) results in 


(E*— 1) 
(9a) 
W 46 logE 


For values of E> 10, (9a) may be simplified to 


E 1 


1+ W 46 log E 

From this equation it is seen that reduction of 
H_ W to zero value gives rise to the original Bowman- 
CICHELLI equation for total reflux. 

In order to make a rapid solution of eq. (9a-9b) 
possible, a number of s values were calculated as a 
function of different log 2 values and H/W ratios. 
The calculated values are given in Table 1, from 
which the graph in Fig. I-A has been constructed. 
In this graph log =A has been related to the 
separating factor S(= 100s) with H W as a para- 
meter. The values of log E=A obtained may be 
converted to the minimum plate number by 


A 
N,, = (10) 


The solution of this equation is simplified by the 
use of the alignment chart in Fig. I-B. 
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(2) Minimum reflux ratio 
For the calculation of the minimum reflux ratio the 
starting point is given again by eq. (4). The relation 


A 


al 
a Been ow 


8 


a 
= 
= 


8 


ITED! 


TT 


TT 


Q? 


Diagrams I-A and I-B for calculation of minimum plate 
number 


between top and bottom composition needed to solve 
this equation is obtained from the condition of 
minimum reflux, ¢.g. infinite number of plates in the 
column (H.E.T.P. infinitely small). According to 
this stipulation the equilibrium line and operating 
line in the McCabe Thiele diagram intersect; and at 


the point of intersection the following relations hold: 


ax 
1+(« 


(Il) 


and 
Ry, 


+ 


X = 


Combination of these equations gives the relations 


between y and x, 
(Ry + 


(I3a) 


and 
dy 


(Ru + 1) 
dx R 


In order to obtain the relation between holdup 
composition and bottom composition, it may be 
pointed out that at the intersection of equilibrium 
line and operating line, the liquid is in equilibrium 
with the vapour ascending from the reboiler. The 
composition at this poigt therefore equals the still 
composition. Now, as the plate number in the column 
is infinite, we may put: 


(I3b) 


dz, 


dx (14) 


For solution of eq. (4), this relation is rewritten: 
H 


d 
B=(y—2) 


B=s(1+ (15a) 


and 
(l5b) 


By substitution of (13a) and (13b), &,, may be 
solved from these equations. However, this method 
gives rise to rather complex relations, which cannot 
be handled easily. Therefore, a numerical solution 
was chosen, which was obtained as follows: 

Eq.(l5a) expresses the magnitude B to be the 
separating factor corrected for holdup. If H W 0, 
B equals s. The relation between B, the quality of 
separation and holdup quotient has been given in 
diagram 

Moreover, magnitude B was calculated with 
eq. (15b), employing a number of «-values and 
minimum reflux ratios. First, 2 was solved from 
eq. (13a) for «-values of 1-05, 1-10, 1-25, 1-50 and 2-00 
and a number of minimum reflux ratios (2-200) in 
the case of y=0-5. 

These x-values were substituted in eq. (13b), 
from which dydx values were obtained. Finally, 
the x and dydz values were used for calculation of 
B according to eq. (15b). 

The data obtained are collected in Table 2. In 
order to construct a general relation between B, R,, 


10 


x 
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and x, reciprocals of s and R,, were plotted on a large- Table 2. Calculated B values in column with infinite plates 
scale graph. The curves for each «-value were drawn when distillate composition y = 0-5 
through the origin at which R,,= 00, s= oo. From 


Ry dy/dzx B 
these curves s-values up to a maximum of 25 were 
read off. 1-05 20 0-216 1-437 0-310 
40 0-314 2-260 0-710 
| 3 2-092 
200 452 10-041 4-539 
I-A | 1-10 1 1-445 
' 20 2-268 0-714 
100 0-453 10-063 4559 
‘ Y, 1-25 5 0-262 1-681 0-440 
10 2-754 0-956 
20 | 0-412 5-132 2-114 
+ 40 0-453 10-066 4-560 
| LY | 1-50 2 0-253 1-565 0-396 
+ + 5 0-353 2-810 0-992 
Q tH 10 0-413 | 5-153 2-128 
= 2-00 2 0-343 2-482 0-851 
5 O-417 5-231 2-181 
0 0-453 10-069 4-561 
43 
x 4 
1§ x,—causing the numerator in eq. (13a) and (13b) to 
= be approximately equal to 1—relation (15) reduces to 
45 
s=y(R,,+1)(«a—1) (15a) 


Vs 


}— 
Diagrams IT-A and II-B for calculation of minimum + Points, interpolated trom a 
reflux ratio lineair scale t/s vs 1/Rm plot 


Construction of a simple alignment chart based 
on the data obtained was not possible as the relations — Som oo on Yam! 
between Is and 1/R,, were not linear. The same Fig. 1. Plot of Is vs 1)R,,. 

was found when plotting log 1/s vs log 1/R,,. However, 

by deforming the log 1/R,, scale, a linear relationship — from which, when y= 0-5: 
between log 1's and /(log 1/R,,) could be obtained. 


2 
The result is shown in Fig.1, from which the R,~, ~: (15b) 
alignment chart in Fig. II-B has been constructed. 
It may be noted that in the case of negligible This is the Bowman-CIcHELLI equation for mini- 


holdup and low value of the bottom composition, mum reflux. It holds good at s-values above about 


11 


20; at lower s-values it tends to give too low reflux 
ratios, which may be checked by the alignment 
chart IL-B. 


(3) Relation between the minimum requirements and 
plate number and reflux ratio to be used. 


In actual practice distillation columns wiht an infinite 
number of plates are not possible as a result of which 
the reflux ratio has to be chosen higher than the 


calculated minimum value. It will b~ clear that 
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followed by averaging by graphical integration. This 
gave the value of dx, dy which, upon multiplication 
by the dy dx value for the total plate number involved, 
furnished the value dx,/dx desired, 

The basic values of (y—a), dydx and dx, dx 
obtained in this way are collected in Table 3. 

In determining the appropriate relations between 
RR, and N_N,, the following procedure was followed: 

Five ratios of holdup/residue were chosen, namely, 
HW-0, 0-05, 0-1, 0-2 and 0-5. For each value 


Table 3. Values of (y— x), dydx and dx, dx, used in solution of eq. (4) 


x 


dx, 


dy 


dx dx 


25 


dy 


dx 


dz 


20 0-155 115 5 
40 0-200 1-31 10 0-245 
80 0216 1-44 20 0-262 
160 0-216 1-44 - 40 0-262 
40 20 0-186 1-18 10 5 0-217 
40 0-268 1-49 10 0-304 
80 0-309 214 20 O-347 
160 O314 2-26 O-347 
100 20 0-209 1-20 20 5 0-230 
40 0-325 1-79 1-0 10 347 
80 0-399 3-7 2-18 20 0-409 
160 0-410 | 2-17 40 0-412 
200 20 218 1-25 — 40 5 0-244 
40 0-347 20 1-60 10 0-373 
80 (437 57 3-06 20 447 
0-452 341 0-453 


1-20 ‘ 
1-40 5 0337 2.06 

1-68 10 248 — 
1-68 - 5 2 | 0281 1-43 1-28 
1-22 5 O401 3-42 2-19 
1-82 10 O417 2-27 
2-75 10 2 0200 1-56 1-38 
2-75 5 0430 472 2-82 
1:30 1-20 10 0453) (101 3-68 
214 1-62 

4-07 2-15 

513 1-89 

131 1-19 

2-38 1-74 

6-37 3-06 

2-92 


with increasing NN, ratio, the ratio RR, will 
approach unity. 

Determination of the relation between R R,, and 
implies again the solution of eq. (4). The 
relations needed for dydx and dx, dx as a function 
of x, however, are so complicated that an analytical 
solution of the problem is not feasible. 

Therefore, the problem was attacked by means 
of graphical calculation methods. 

As a starting point y—~- relations were calculated 
with a modified equation according to Smoker [16]. 
a-values of 1-05, 1-25 and 2-0 were used and a large 
number of plate values and reflux ratios. The y—x 
relations obtained were plotted and at y—0-5 the 
values of and dy dx were read off. 

In addition, the values of dx,/dx had to be cal- 
culated. These were obtained by plotting the dx dy 
values vs the plate number at a fixed reflux ratio, 


substitution of the values from Table 3 in eq. (4) led 
to a series of values of s as a function of R and N. 
Two graphs were made from these relations. In the 
first 1s was plotted vs 1N with R as parameter, 
Drawing the curves obtained to the limiting value 
of Ls with 1. N —0 enabled reading off of 1/s-values 
for every desired plate number. 

In the second plot 1s was plotted rs 1 R, N now 
being the parameter, using data from the first plot 
too. Drawing the curves obtained to the limiting 
1/s-values at 1/R=0 made readings at high reflux 
ratios possible. 

From this graph data of 1/s-values (and con- 
sequently, s-values) were taken at a series of pre- 
determined reflux ratios, involving a series of plate 
numbers, These data were finally replotted in a log S 
vs log N graph, R being the parameter. Readings of 
plate number at «= 20, 10, 5, 1 and 0-5 were taken 
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Table 4. (Continued) 
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at the predetermined reflux ratios mentioned above. 
In this way for each s-value a series of plate numbers 
and reflux ratios were obtained which formed the 
starting point for the final correlation. The data are 
collected in Table 4 together with the minimum 
requirements R,, and N,,. 

In studying these data it was found that they 
could be represented by the following relation between 
R, R,,, N and N,: 

R 


R—R, |N 


(16) 


Relation between power coefficient n 
and constant C in eq. (16). 


Fig. 2. 


This equation fulfils the boundary conditions, 
namely at R= R,,, N= oo and at N=N,, R= ow. 
The power coefficient n was found only to be a 
function of C, which relation is shown in Fig. 2. 


Again the constant C is a function of the pole height 


s and the holdup quotient HW. No injluence of 
a-value was found. 


14 


The relation between constant C and the separat- 
ing factor S = 1008 and is shown in Fig. III-A. 
Eq. (16) was modified in an alignment chart, given in 
Fig. III-B. 

In order to give an impression of the accuracy of 
the correlation eq.(16), the data in Table 4 were 
recalculated with the help of Fig. III-A and III-B. 
On the one hand, using the reflux ratios and R,, and 
N, as the starting point, N was calculated; on the 
other hand, using the V-values as the starting point, 
R was calculated. The data obtained are shown in 
Fig. 3 as a function of the N and R values given. 
Agreement between calculated and given N values is 
generally quite good. The largest discrepancies are 
found with the calculated reflux ratios. This is due 
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to the circumstance that the N/N,, ratios encountered 
are generally small, resulting in a large variation in 
R/R, when N/N, varies only by a small percentage. 


a 


Te 


Diagrams III-A and III-B for calculation of actual reflux 
and plate number. 


This may be checked directly with Fig. III-B. It 
is felt, however, that the accuracy of the calculation 
method is within about 10%. 


APPLICATION OF THE RELATIONS OBTAINED 


Calculation of the minimum plate and reflux require- 
ments besides the actual plate number and reflux 
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ratio to be applied may be rapidly made by using the 
diagrams in Fig. I-III. 

The following has to be known: 

1. Relative volatility x of the components to be 
separated. 

2. The sharpness of separation desired. The slope 
of the distillation curve to be obtained is related to 
the purpose of the distillation. Expressing this slope 
as mol.-% volatile component in distillate per mol.-°% 


R values table 
1000 


Fig. 3. Test of eq. (16) (diagrams III-A and III-B). 


residue the following values may be used as an 
orientation : 
Mode of separation 
sharp 


normal 


bad 


Slope 

20-25 
8-12* 
3-5 


3. The residue W and holdup, H. Where the 
distillation is carried out with appropriate sharpness, 
the sum of holdup and residue, when distillate com- 
position is 50°, may be put approximately equal to 
the amount of heavy component in the original feed. 
Expressing holdup and content of heavy component 
as mol.-% of the original charge the residue is ob- 
tained as the difference between these two quantities. 

The values of S and H/W can now be calculated. 

Using the values obtained, the constants A, B and 
Care read off in Fig.I1-A, and III-A, respec- 
tively. The constants A, B and C are then used in 
combination with the relative volatility « in the 
alignment charts I-B, II-B and III-B to obtain N,,, 
R,, and the proper N/N,, or R R,, ratio. 

* Nearly equivalent to the standard separation of binary 
mixtures containing 50% volatile component according to 


Rose {10}. In this separation the first 40% of cistillate 
obtained has an average purity of about 95%. 
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The method described has been derived for binary 
mixtures only. However, in a distillation of a mixture 
with more than two components when the components 
are obtained ir reasonable purity (a mid fraction not 
being contaminated both with adjacent more volatile 
and heavier components) the plate and reflux require- 
ments may likewise approximately be found by the 
given method. 

To this end, the multicomponent system is divided 
into a number of binary systems, in which the residue 


55 ml under the conditions applied. The following 
mixtures were distilled under various conditions: 
methyleyclopentane-cyclohexane (mep-ch) 
n-heptane-methyleyclohexane (nh—mch) 
meth yleyelohexane-dimeth yleyelohexane*(mch-dmch) 
cyclohexane- meth yleyclohexane-dimethyleyclohexane 
(ch-meh-dmch) 

The values of dydw, H and W were determined 
from the distillation curves obtained; the data are 
given in Table 5. This table includes values of R,,, 


Table 5, Experimental verification of the calculation method 


Slope 
Composition -% 
moal.-% on. 


24mep; 28 
75 nh; 25 meh 1s 
30 ch; 30 mch; (a) 
40 dmch 14(b) 
35 ch; 30 mch; 23(a) 
35 dmch 22 (b) 


29 mech; 7idmch 46 5 0-53 1- 


0-075 1- 


N 
exp. 


14 12% 
11-5t 

97 

9-3 

| 

87 


9-2 


VOL. 


1952 


35 ch; 30 mch; 18(a) . 2010 18 . 7-2 
35 dmch 17 (b) 0-075 1-7 
75 nh; 25 mch 3-0 . 57 O35 1-07 O08 27 1-75 1-65 
2800 26mep; 74ch 11 0-085 1-3 8-4 215145 51 35 28 1-15 


+ calculated with eq. (9a); (a) at separation ch-mch; (6) at separation mch—dmch. 


W is again found as the difference between content 
of heavier component of the binary system considered 
(expressed in mol.-°%o of original charge) and holdup H. 


EXPERIMENTAL VERIFICATION 
or THE CALCULATION METHOD 


A number of batch distillation experiments with 
binary and ternary systems were carried out in order 
to obtain an experimental verification of the calcu- 
lation method outlined above. These distillations were 
carried out in Oldershaw sieve-plate columns consist- 
ing of one section (30 actual plates) and three sections 
(90 actual plates), respectively. Reflux ratio control 
was of the intermittent vapour dividing type. The 
apparatus has been extensively described by CoLLins 
and Lantz [5}. 

Calibrations with the normal heptane methyl- 
cyclohexane binary system showed the stationary 
theoretical plate number of the one-section column 
and three-section column to be about 17 and 51, 
respectively. Holdup per section amounted to about 


N,,, and reflux ratio to be used, calculated with the 
diagrams I-III. Relative volatilities x were taken 
from vapour pressure data from A.P.I. project 44 of 
National Bureau of Standards, Washington. 

The experimental reflux ratios in Table 5 under 
heading I show a reasonable agreement with the 
calculated ones. This is also the case with the separa- 
tions in the ternary system ch-mch-dmeh, which 
proves that the calculation method may be applied 
to this type of systems, too. 

The conclusion is reached that in these experiments 
the conditions of equilibrium between the contents of 
the column and reboiler was met. 

However, another picture is obtained from the 
experiments under heading II in Table 5. The reflux 
ratios applied are in all cases much lower than the 
calculated values. With the separation of mep-ch in 
the column with 51 theoretical plates the reflux ratio 
applied is even appreciably lower than the minimum 
value. 


* Mixture of mainly trans 1-3 dmch and trans 1-4 dmch 


Exp. Cha n N , | 
Group “ye mol.-% mol-% A ix, ( | cate.| cap. 
1 930 59 69-1 0-085 13 190 2-05 16 17 1:35 27 21 32 30 
2 930 58 72 080 107 #13 0-23 11 7615 11 27 
3 730 76 224 035 18 380 51 17 175 15 20 3 
6-5 33-5 020 17 470 57 17-5 17 +18 195 165 29 #£=28 

4 2000 28 272010 18 630 69 18 17 22 26 1313 | 24 28 po 
2-4 326 7:7 23 17 «1-95 30 4135 31 £28 
6 9 #44 
Il 23 
7 38) O30 
8 67 30 

16 
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The results of these experiments suggest that the 
concentrations in the column were not in equilibrium 
with the reboiler concentration. This effects a better 
separation than may be expected when equilibrium 
exists. The beneficial effect of holdup in batch 
distillations was also observed by CoL_BuRN and 
Stearns [4] and Rose and co-workers [{11, 15]. Rosr, 
Wituiams and Prevost found that at low reflux 
ratios holdup had a beneficial effect on the separation, 
while at high reflux ratios the holdup causes a decrease 
in sharpness of separation. The critical reflux ratio 
at which holdup had no influence upon separation 
quality increased with the plate number of the 
distillation column used. 


In our case the deviation from equilibrium in the 
column may be due to 


(a) the increased rate of distillation, caused by 
the lower reflux ratio in exp. 5 and 6. This results 
in a higher rate of change in concentrations in the 


column. 


(6) the increased theoretical plate number in exp. 7 
and 8, causing a slow rate of equilibration of the 
column. Calculation of the time needed to reack a 
certain fraction of the stationary plate number accord- 
ing to the method of Bera and James [1] shows for 
instance an equilibration time 20 times larger for the 


5l-plate column than for the 17-plate column. 


The experimental investigation leads to the con- 
clusion that under certain conditions the calculated 
reflux and plate requirements agree with the values 
to be applied in practice. In some cases, however, 
due to insufficient equilibrium conditions in’ the 
column, the values to be used in practice are lower 
than those theoretically predicted, 
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NOTATION 
A = constant in eq.(10), A = log 
B — constant in eq. (15), B=s(l + 


C — constant in eq. (16) 

E fractionating factor, 

H - holdup of fractionating column, expressed as % of 
original charge; mol.-% 


N — theoretical plate number 


N minimum number theoretical plates 


m 
n = power coefficient in eq. (16) 
R = reflux ratio, quotient of reflux flow and product flow 


R,, ~ minimum reflux ratio 
dy 
aw 


mol. fraction 


SN — separating factor W 100; mol.-% 
a 


pole height w 

dw 

W — residue in reboiler, expressed as % of original charge; 
mol.-°%o 


X — composition of vapour, in equilibrium with reboiler 
liquid (volatile component); mol. fraction 


x — composition of reboiler liquid (volatile component); 
mol. fraction 


x), ~ average composition of column holdup (volatile com- 


ponent); mol. fraction 


y ~ composition of distillate (volatile component); 
mol. fraction 


relative volatility 
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The statistical interpretation of thermodynamic concepts 
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Summary—The spread of molecular velocities, which was used by MAxwett and Bourzmann for the 
kinetic theory, is shown to provide a simple pictorial representation of thermodynamic quantities. 


Résumé— L’ auteur montre que la répartition des vitesses moléculaires, déja utilisée par MAXWELL et Bourz- 
MANN dans la théorique cinétique des gaz, peut fournir une représentation imagée des grandeurs thermodyna- 


miques. 


The teaching of thermodynamics to chemical engineer- 
ing students is particularly difficult because there 
are two classical approaches to the subject. The 
one is in the design and operation of heat engines 
and the other is in the physical chemistry of reaction 
processes, 

In many cases these approaches follow the same 
argument, although often couched in different 
language and using different symbols. Considerable 
strides have been made to overcome these purely 
mechanical difficulties by the trend towards standardi- 
sation of symbols and the publishing of books in 
which the subject is treated from both points of view. 

The fundamental difficulty still remains however 
that the student has to acknowledge the fact of the 
second law of thermodynamics without any pictorial 
or other more satisfactory explanation. The second 
law can be, and has been, stated in a legion of different 
ways but for the engineer it is best expressed in the 
simple form that out of heat energy available (such 
as obtained in the burning of coal) only a certain 
proportion can be obtained as work energy (such 
as electrical or mechanical power) and that this 
proportion depends on the upper and lower tem- 
peratures of the particular cycle used (that is the 
steam and condenser water temperatures). 

The concept of entropy is also far from clear. 
He has to accept it as a mathematical entity without 
an understandable physical background. When 
R. CLausivs wrote his paper on “ Different forms of 
the principal equation of the mechanical heat theory” 
{1} he realised that the quantity [dQ/T (where dQ 
is the heat absorbed by a system at an absolute 
temperature 7’) was worthy of a special designation 
since it was an exact differential depending on the 
state of a substance and not on the method by which 
it reached that state. He called this quantity “S”’ 
and said: 

,, Sucht man fiir S einen bezeichnenden Namen 
so kénnte man, ahnlich wie von der GréBe U gesagt 


ist, sie sei der Wirme- und Wiirmeinhalt des Kérpers, 
von der Grébe S sagen, sie sei der Verwandlungs- 
inhalt des Kérpers.* 

Verwandlung” is translatable into English as 
“alteration, transformation, 
metamorphosis’ “changing one 
figure, or quantity, into another equal to it.’’ CLav- 
situs selected the Greek 7 room) and coined the word 
“Entropie”’. 

It is perhaps idle to examine this definition to 
try to discover any more in it than is contained in 
the mathematics. It is perhaps too much to expect 
that a mathematician is a better exponent of language 
than mathematics. 

When J. Wittarp Gripes wrote his paper ‘On 
the Equilibrium of Heterogeneous Substances” [2] 
in 1876, he considered, among other subjects ‘‘Con- 
siderations relating to the increase of Entropy due 
to the Mixture of gases by Diffusion.’’ He then made 
the following statement “The impossibility of an 
uncompensated decrease in entropy seems to be 
reduced to improbability.”” This sentence is by no 
means a model of lucid English but from the context 
seems to mean that ‘To separate a mixture of two 
gases without external influence seems improbable 
but not impossible.” The sentence is important 
however since it gave BourzMann the suggestion 
of relating entropy to probability. He quotes it, for 
example, in the heading of the second part of his 
book on ,,Vorlesungen iiber Gastheorie’ in 1896 
where he discusses probability under the section 
,,Mathematische Bedeutung der GréBe 

Let us now for a moment turn our attention to 
the kinetic theory. This theory has developed our 
knowledge of the behaviour of molecules on a physi- 
cal basis but has not, so far, been utilised to give us 
an adequate explanation of these thermodynamic 
problems. At its inception the statistical distribution 
of molecules was assumed although it is only com- 
paratively recently, by the researches of O. Stern, 


change, conversion, 


or in mathematics as 
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I. F. ZartMan and others, that it has been shown 
experimentally that the distribution of molecular 
velocities does indeed follow a statistical law. 

The marriage of these two methods of examining 
the state of matter is eminently desirable and it is 
the aim of this paper to show, without recourse to the 
necessary background mathematics, how this can 
be done [3]. 

In plotting the distribution of molecular velo- 
cities from the experimentally determined trace of 
bismuth molecules on the glass support the ordinate 
is obviously the number of molecules with a given 
velocity. This is not the only way of plotting a 
statistical distribution since the surface area of an 
emulsion, or the weight distribution of a powder, 
can be derived from the count of the number of 
particles with any given diameter d and multiplied by 
d? in the one case and d* in the other. Similarly instead 
of plotting number » against a function of the velo- 
city, it is possible to plot nv*® which gives a measure 
of the kinetic energy of the molecular species with a 
given molecular velocity. Thus the spread of mole- 
cular energy with molecular velocity can be represen- 
ted by the statistical curves shown in Fig. 1. It is 
best to use the term a “function of the velocity” 
for the abscissae since it will then take care of the 
various rotational and vibrational motions charac- 
teristic of the liquid and solid states of matter. 

Now a characteristic of a statistical distribution 
is that it requires defining by both the mean and 
the spread and it can be shown by MAXxwetts cal- 
culations that the energy is greater as the spread 
increases. We have therefore a method of representing 
entropy as a mathematical entity connected with 
statistics and which, for a given temperature, can 
have an infinite number of possible values. This con- 
cept is illustrated in Fig.1 where the system has the 
same mean but different entropies or spreads. In 
passing it might be mentioned that in calculating the 
entropy of a solid the only physical measurements 
used are those of specific heat and temperature; 
for a liquid the latent heat of melting has to be in- 
cluded; and for a gas the latent heat of vaporisation. 
The generalisation known as Trovutons rule that the 
entropy of vaporisation is approximately constant 
on this concept means that the spread of the distri- 
bution on going from a liquid toa gas is increased ap- 
proximately the same per molecule for all substances. 

The statistical interpretation gives a clear idea 
of the reason for the second law of thermodynamics. 
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The diagram for two different temperatures is shown 
in Fig. 2 where it will be seen that the two curves 
overlap. The molecules with the energy content 
represented by the hatched area are common to 
both temperatures and therefore do not give up 
energy in changing from the higher temperature to 
the lower. All the molecules in the unhatched area of 
B are available to giving up part of their energy 


Sw) 
Fig. 1. 

to those in the unhatched area of A which are lower 
on the velocity scale. As the temperatures draw 
apart so do the curves A and B and the hatched 
area becomes smaller. Conversely when the curves 
A and B draw together the hatched area becomes 
larger and less of the total energy of B can be con- 
verted into work energy when the system changes 
to the temperature state indicated by A. 


nv? 


dh 


Fig. 2. 


When chemical reactions are considered there are 
two main quantities, the Free Energy and the Heat 
of the Reaction, to be explained. The term “Free 
Energy”’ comes from consideration of 
“Bound Energy” and it is curious that the latter 
word has been discarded but not the former, which 
has been derived from it. The equation relating these 
two quantities of Free Energy and Heat of Reaction is: 


AG=AH—TAS. 


For a reaction to take place spontaneously there 
must be a decrease of Free Energy of the system but 
the reaction can be either “exothermic” or “‘endo- 
thermic.” From this it follows that the entropy S 
may either decrease or increase and this can be 
represented on the statistical theory by a change in 
the spread between the distributions representing 
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the reactors and the reactants. If the entropy of 
the reactants is sufficiently smaller than that of the 
reactors the sign for the Heat of the Reaction will be 
the opposite to that of the Free Energy. 

One of the main subjects of heat engine theory 
is concerned with the adiabatic or is entropic expan- 
sion of gases. According to the kinetic theory there 
is an identity between the mean velocity U’ of the 
system and its pressure which is given by: 


PV=RT=4MNU?. 
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fore a number of degrees of freedom allowed to a gas 
under these conditions without altering the spread of 
the distribution. The fact that the velocities are 
temporarily unidirectional instead of random would 
indicate that they are not available to exert pressure 
in the form in which it is customarily measured, 
The treatment of this theme is, of necessity, 
brief and incomplete but it is felt that the picture 
here presented is a sufficiently clear one to be of 
considerable value in explaining the problems of 
thermodynamics in a more readily assimilated form. 


The mean velocity /’ is of the same order as the 
velocity of sound in the medium and, according to 
adiabatic expansion theory, is the limiting velocity 
that a gas can attain when allowed to expand adia- 
batically through a simple nozzle. There are there- 
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Summary —Formulae are derived, concerning the separation of solids from liquid by a continuous filtration 
procedure, consisting of the filtration proper and the washing of the filter cake by a suitable wash solvent. 

In the derivation of the formulae it has been assumed that the mixture of solids and liquid to be filtered 
is prediluted with a low-viscosity liquid in order to make possible or to speed up filtration. However, 
the formulae can easily be adapted to cases where no predilution is used. 

The effect of the process variables on cake impurity and liquid yield is shown. 
examples have been tabulated and plotted for practical purposes. 

The beneficial effect of filtrate recirculation on the degree of separation, already known in filtration 
practice, has been worked out quantitatively in formulae, tables and graphs. 

On application of the formulae to the dewaxing of lubricating oils with the aid of butanone-benzole 
mixtures, used both as diluent and as wash solvent, good agreement with practical results was obtained, 


Some numerical 


Résumé—L’auteur établit des formules relatives & la separation liquide-solide dans le cas d'une opération 
continue, comportant & la fois la filtration proprement dite et le lavage du tourteau par un solvant con- 
venable. 

Pour lestablissement des ces formules, on suppose que le mélange initial liquide-solide est dilué au 
préalable avec un liquide de faible viscosité pour accélérer ou méme rendre possible la filtration. Toutefois, 
les formules peuvent étre généralisées aisement dans labsence du toute dilution préalable. 

Des exemples numériques, présentés dans forme de formules, tableaux et graphiques, montrent 
leffet des divers facteurs de lopération deu la teneur en impuretés du tourteau et sur le rendement en 
liquide. 

L’auteur etudie quantitativement (formules, tableaux, graphiques) le recyclage du filtrat, dont leffet 
avantageux sur le degré de séparation est déja étalli par la pratique industrielle. 

L’auteur presente des resultats expérimentaux en bon accord avec des formules dans le cas du déparaf- 
finage des huiles de grainage au moyen de mélanges benzé 1e-butanone, utilisés 4 la fois comme deluants et 


comme liquides de lavage. 


1. Iwrropvucrion (a) The liquid, which is retained in the, generally 
When the separation of a solid and a liquid phase highly porous, filter cake, can be separated from the 
is desired [subsequent to filtration (or centrifuging)|, solids therein by distillation or drying, owing to a suffi- 
two cases may be distinguished: cient difference in volatility between solid and liquid. 
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(b) Separation of solids and liquid cannot be 
effected by distillation or drying, e.g. because both 
phases comprise components of the same volatility 
such as in the de-waxing of lubricating oils. 

In the latter case the liquid, retained in the pores 
of the cake by capillary forces, is partly or wholly 
removed by a wash liquid, which is so chosen that it 
can, in turn, easily be separated from the washed 
cake and the wash filtrate, e.g. by distillation. 

Addition of a low-viscosity liquid to the suspension 
to be filtered in order to make possible or to speed up 
filtration, is often necessary or desirable. 

The following deals with the most general case that 

1. a low-viscosity liquid, hereafter called solvent 
or primary solvent, is added to the suspension to be 
filtered, and 

2. a wash liquid is applied to the filter cake. 
Special cases (no primary solvent and/or no wash 
liquid applied) can, however, easily be derived from the 
general treatment of the problem. The primary 
solvent is so chosen that at the filtration temperature 
it mixes completely with the original liquid and dis- 
solves only to a minor extent the solids present in 
the suspension. Besides, it must be easily recoverable 
from filtrate and filter cake. 

As the wash solvent has to meet the same require- 
ments, in most cases the same solvent is used for 
diluting the original suspension and for washing. 
Ideal washing (pure displacement of liquid in the 
pores by wash liquid) would only require an amount 
of wash liquid equal to the pore volume of the cake 
to replace the liquid in the cake completely by wash 
liquid. Although ideal washing is theoretically not 
possible and hence never met in practice, yet the 
volume ratio of wash liquid to liquid to be washed 
out (= pore volume) determines washing efficiency. 

Calculation of the degree of separation between 
solids and liquid in continuous filtration therefore 
requires an expression for the relation between the 
volume ratio of wash liquid to pore volume and the 
proportion of liquid washed out. Besides, the amount 
of wash liquid that can be applied to the filter cake 
during the available wash time, based on the pore 
volume of the cake, must be known. 

If it is assumed that the primary filtrate has 
the same composition as the liquid present in the cake 
before washing, and that the pore volume remains 
constant during washing, it is, in principle, possible 
to derive equations for cake impurity and the yield 
of original liquid as a function of the process variables. 
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2. DeERIVATION OF THE Basic EQuaTIon 


The following equations, although discussed mainly 
in connection with rotary drum filtration, are applic- 
able to any other kind of continuous filtration, as 
defined in Section 2 of our previous paper [1]. 

To each segment, with surface area 4A, of a 
rotary drum filter the batch filtration equation for 
constant pressure 


dq p.4AA 
dt kq 


(1) 


+r] 


is applied for the filtration proper (see Section 1 of 
the previous paper). 
Integration of (1) leads to: 
1 kg? : 
+nHrq=pT.AA, (2) 
where 7’ is the filtration time of the segment per 
revolution of the filter. (2) can be transformed into 
an expression for the mean filtration rate: 


nls +r) 


(3) 


The difference between (3) and (1) is only that for the 
mean filtration rate the sum of half the cake resist- 
ance at the end of filtration and the cloth resistance 
must be inserted, whereas for the momentary filtration 
rate it is the sum of momentary cake resistance and 
cloth resistance. 

Although the supposition on which (1) is based, 
viz. that cake resistance increases linearly with the 
amount of filtrate obtained, has proved its validity in 
practice, it is worthwhile to consider the effect of 
deviations from this supposition on the mean filtration 
rate. 

If the batch filtration satisfies the equation 


dq p.AA 
dt (la) 
ad 
the mean filtration rate is found to be: 
£p.4A 
(3a) 


1 
iat!) 


As n usually will not differ very much from the value 


n= 1, the factor 


the value 4, so that also for deviations from (1) it is 
reasonably accurate to say that the mean filtration 
rate is determined by the sum of half the cake resist- 
ance at the end of filtration and the cloth resistance. 


will not differ very much from 
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After filtration the cake is washed during the time 
interval 7. The resistance to be overcome by the 


wash liquid is constant and equal to| a, ae r), so that 
the mean wash rate is expressed by: 


(4) 
Te kq +r 
"(44 


if the pressure difference, p, is the same during washing 
and filtration. 

As the viscosity of the wash liquid generally 
differs from that of the liquid to be washed out, the 
flow rate during washing is not necessarily constant. 
Therefore 7* means a viscosity intermediate between 
7 and »,, the viscosities of filtrate and wash liquid, 
respectively. 

For the surface area 1A and hence for the total 
filter area A, the ratio of wash liquid to primary 
filtrate sucked off during a certain time interval is 
found by dividing (4) by (2), giving: 

. 
qu 2 A (5) 
&,, 
A 
[qe, the amount of wash liquid applied, is equal to 
the amount of wash filtrate obtained, if, as supposed, 
the pore volume of the cake remains constant during 
washing. | 
kq 
Ar 
cloth resistance at the moment that washing is started) 
(5) becomes: 


By introduction of «= (ratio of cake and 


ja+l 
a+l° 


de Te 


(6) 


From this equation it follows that q,/q is twice as 
high for «=0 as for «oe, in accordance with the 
results of the preceding article, in which it has been 
shown, that by increasing the speed of a rotary drum 
filter, i.e. by decreasing x, at most twice as much 
wash liquid passes through the same amount of cake. 
Te n jat+l. 
T a+1 
put — x, and both terms are divided by g,, the amount 
of liquid present in the cake after filtration (— pore 
volume of the cake), (6) can be written as: 


When, for reasons of simplicity, 


(7) 


Thus an expression has been found for the ratio of 
wash liquid to liquid to be washed out in continuous 
filtration. q q,, the ratio of filtrate to liquid to be washed 
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out, is equal to the ratio of original liquid in the 
filtrate to original liquid in the cake before washing, 
because, if no preferential adsorption takes place, 
both liquids have the same composition. 

When the fraction of original liquid removed in 
the primary filtrate is designated by O, this means 
that 


q oO 
oO’ (8) 
so that the basic wash equation is: 
qu _ O 
"1-0" (9) 


For the separation of a suspension into solids and 
liquid on a rotary drum filter it appears that if x is 
considered as a constant, q¢,/q, is a simple function 
of O only. O itself, however, is a function of the 
amount of primary solvent used, solid content of 
the mixture to be separated and porosity of the cake; 
the structure of this function will be derived in the 
next section. 


3. Tue Bastc Equation ty Terms or SLURRY 
AND CAKE VARIABLES 

For the purpose of working out eq. (9) we introduce 

the following quantities: 


F — solid content of base material to be separated 

S=ratio of dilution solvent to base material 
(= solvent ratio) 

e = liquid content of cake (= porosity) 


These three quantities, as well as the quantities that 
will be introduced later, are expressed in volumes. 

By base material is always meant the solid-liquid 
mixture not yet diluted with solvent. The diluted 
mixture is normally called slurry. 

The cake before washing consists of: 

F solids and," liquid, provided the solubility 
of the solids in the solvent can be neglected. 


The liquid in the cake can be divided into: 
eF i1—F F 
ei—-F+8 
solvent. 

The total amount of original liquid present in the 
base material equals (1—F), so that (l1—O), the 
fraction of original liquid left behind in the cake after 
filtration, but before washing, is given by: 


original liquid and 


eF 1 


l—e 1-F+8 


(10) 


Vol. I 
No. 1 — 1951 


and 


by: 


0 
1—O eF 


(9) can now be written as: 


(l—e)(1+S8)—F 
x 
eF 


(12) 
By introduction of S’ = 
liquid to solids in the slurry, this equation is trans- 
formed into: 


, the ratio of 


qu l—e 
= 2/( S’—1). 


(13) 


In this way q,/q. has been made independent of F. 

From equations (12) and (13) several conclusions 
can be drawn: 

1. The influence of the cake porosity ¢ on q,/q., 
and hence on the amount of original liquid removed 
from the cake by washing, is extremely high because 
of the presence of the factor 

This means that it is important, when the solid 
particles in the slurry have been obtained by crystal- 
lization, so to control crystallization conditions that 
a low porosity of the cake results, if possible without 
simultaneously decreasing cake permeability and thus 
the filtration rate. 

2. A higher primary solvent ratio S increases 
Gul%, for constant values of ¢, F and x. 

3. A higher solid content of the base material F 
lowers q,,/q¢-, When S is kept constant, and has not 
any effect on q,/g., when S’ is kept constant. 
Tw  ja+1 
T a+1 
order to obtain higher values for q,/q, can only be 
effected by increasing the wash time in respect of the 
filtration time. When the same substance is used as 
diluent for filtration and as wash solvent variations 


in »/7* are usually small. 
ja+l 
a+l 
tion of the speed of a rotary drum filter or of the 


cycle time in another continuous filtration operation. 


4. Important variations in x= in 


undergoes relatively small changes on altera- 


4. Cake Impuriry YIELD or OriGinaL Liguip 


From the composition of the slurry that of the cake 
before washing can be calculated, when the cake 
porosity is known. The amount of wash solvent 
that can be applied based on volume of pores, is 
derived from (12). The fraction of the liquid removed 
from the cake by washing, can be designated by 
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f(e/q-) and is equal to the fraction of original liquid 
removed from the cake by washing. 

The composition of the cake after washing can 
now be calculated and from this: 

1. the cake impurity, defined by the percentage of 
original liquid on total amount of solvent-free cake, and 

2. the yield of original liquid, defined by the 
percentage of original liquid, obtained in primary 
filtrate and wash filtrate together. 


Fig. 1. Fraction of liquid washed out f(g, V8- %p Ge- 


In Table 1 a scheme has been given of the express- 
ions required to calculate a continuous filtration. 
Cake impurities and yield of original liquid are given 

J]. expressed in S, the primary solvent ratio, 

2. expressed in P, the total solvent ratio used 
(primary solvent +-wash solvent). 


5. Numerical EXAMPLES 

Tp 
T 
This value may 


For quantitative calculations a value of « = 

» #1 
“n* a+l 4 
occur in rotary drum filtration, if e.g. 7,,/7 = 4 (180° 
submergence of the filter in the slurry, 90° wash 
section), 7 = 7* and % 


will be chosen. 


For /(q»/q-) & curve is used as plotted in Fig. 1. 
The curve is based on the fact that until the moment 
at which an amount of wash filtrate equal to about 
half the volume of pores has been obtained, the wash 
filtrate has the same composition as the primary 
filtrate (pure displacement). Then wash liquid 
“breaks through’’ and the curve deviates from the 
straight line representing pure displacement. This 
curve holds for the washing of cakes obtained in the 
dewaxing process of lubricating oils, where rotary 
drum filters are used and butanone-benzole mixtures 


(1) 
| | | | 
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both as dilution and as wash solvent are applied, at 
least, when cracks in the cake have been avoided. 


Fig.2. Cake impurity before and after washing vs. F for 
constant values of ¢. Total solvent used: 3:5 (see Table 2). 


os; 


Fig. 3. Cake impurity before and after washing vs. F for 
constant values of ¢. Total solvent used: 5-0 (see Table 3). 


Probably, however, the curve has a much more 
general applicability. 

When inserting x=} and the values for 4.) 
from Fig. 1 in the equations tabulated in Table 1, 
results are obtained as listed in Tables 2, 3 and 4. 
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In these tables the influence of cake porosity, solid 
content of the base material, and total solvent used 
can be seen. In Figs. 2, 3, 4 and 5 the results have 
been plotted. 


° 


Fig. 4. Cake impurity before and after washing vs. F for 
constant values of ¢. Total solvent used: 10-0 (see Table 4). 


Fig. 5. Total yield of original liquid vs. F for constant 
values of P and ¢ (see Table 2, 3 and 4). 


Conclusions drawn from the figures are: 

(A) At constant values of porosity and total 
amount of solvent used, the cake impurity before 
washing is lower at higher values of F, the solid 
content of the base material. 
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(B) At constant values of porosity and total 
amount of solvent used, q,/q¢, and hence /(q,/q.) 
decreases with increasing F. 

(C) Maxima in cake impurity after washing may 
oceur as a consequence of (A) and (B) at a certain 
value of F. 

(D) The influence of cake porosity on cake impu- 
rity is so high that the amount of solvent required 
to obtain a cake of prescribed quality, starting from 
a base material with known solid content, cannot be 
predicted without knowing the cake porosity very 
accurately. The same applies to the yield of original 
liquid. 

6. AppLication oF Fivrrate Recirculation 
In filtration practice with rotary drum filters recourse 
is often had to recycling part of the filtrate to the 
filter slurry. By this procedure cakes of higher 
quality are obtained, but at the expense of throughput 
of base material per filter. 

The beneficial effect of filtrate recirculation on 
cake quality can be easily explained for the case where 
the primary solvent ratio, S, is kept constant, for: 

(a) The ratio of solvent to original liquid in the 
cake before washing is the same or higher, in the case 
where a mixture of primary and wash filtrate is 
recirculated. The cake quality before washing is 
therefore at least as good as without recirculation. 

(b) By filtrate recirculation the solid content of 
the slurry is lowered. Hence q,/¢, becomes higher, 
as can already be seen from (12). This means that 
a higher precentage of the original liquid will be 
washed out. 

By keeping the primary solvent ratio constant, as 
assumed above, it is found that the total amount of 
solvent used increases with higher recirculation 
ratios. 

From a practical point of view it is more interest - 
ing to consider filtrate recirculation under the assump- 
tion that the total amount of solvent used is kept 
constant. This means that at higher recirculation 
ratios, the primary solvent ratio decreases. So there 
is a shift from primary solvent to wash solvent, 
which shift can be expected to increase cake quality 
and yield of original liquid. 

Table 5 has been worked out, in analogy to 
Table 1, in order to be able to make numerical 
calculations on the influence of filtrate recirculation. 

R, the amount of recirculation liquid on intake 
of base material, consists of solvent, R,, and ori- 
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ginal liquid, Ry. Ry is computed by the following 
method : 

When R <q, only wash filtrate is recirculated, 
but when R>q, the total amount of wash filtrate is 
recirculated, supplemented by the required amount 
of primary filtrate. 


Example: 
Data: F = 0-20, solid content of base material 
P-—5, total amount of solvent used 
¢ ~ 0-80, porosity of cake 
! constant in basic wash equation 


4? 
= i, recirculation ratio 


z= 

R 
From the equation: 

P= S+2(1+S8+R- i -| in Tabie 5 it follows 
that: 

S=38 

q. = P—S=5-3-8 = 1-2 

qe (l—e)(1+8+R)—F 

qe eF 
From Fig. 1: /(q,,/q,) = 0-90. 


As R=1 and q, =1-2 the recirculating filtrate can be 
supplied entirely by the wash filtrate. 
R,, original liquid in recirculated filtrate, is then found by: 
1 eF i1—F+8, 
1.2 l—e 1—F+S+R 
—> R, = 0-096. 


= 150. 


Cake impurities and total yield of original liquid can now be 
calculated. Results are: 

Cake impurity before washing: 39-1 °% 

Cake impurity after washing: 6-0% 

Total yield of original liquid: 98-4% 


Compared with the figures in Table 3 it can be seen that 
filtrate recirculation has lowered the cake impurity after 
washing from 9-2 to 6-0% and has increased the yield from 
97-5 to 98-4%. 


Other numerical calculations have been tabulated 
in Tables 6, 7 and 8 and plotted in Figs. 6, 7, 8 
and 9. 


7. CoNcLUDING REMARKS 


The applicability of the formulae given in this article 
may be doubted, owing to the yreat number of 
assumptions required to reach the ultimate aim, viz. 
computation of liquid yields and cake impurities in 
continuous filtration. However, some practical expe- 
rience with the formulae was gained in applying 
them to the process of solvent dewaxing lubricating 
oils. The process was carried out with butanone- 
benzole mixtures both as dilution and as wash solvent. 
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Fig. 6. Cake impurity vs. F for constant values of R 
(see Table 6). ¢ = 0-80; P=3-5; x=}. 


%o 


2s 


os 

F 

Fig. 7. Cake impurity vs. F for constant values of R 
(see Table 7). ¢ = 0-80; P=5; r=}. 


For the separation of the wax from the slurry Dorr- 
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(a) The constant x in the basic wash equation is 
lower than would be expected from the ratio of wash 


Cake impurity 
Yo 


= 


Fig. 8. Cake impurity vs. F for constant values of R 
(see Table 7 and 8). ¢ = 0-80; P= 5; a=} and ¢= 0-90; 
P =5; #= }j. 


—— €2080 
€=0.60 
€20.90 P=5 


| | 
Fig. 9. Total yield of original liquid vs. 
values of P, ¢ and R (see Table 6, 7 and 8). 


F 
F for constant 


section to filtration section. This may be due to the 


OLIVER rotary drum filters were used. It was found fact that the wax cake is washed by spraying solvent 


that: 
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Table 8. Influence of filtrate recirculation on cake impurity and yield of original liquid for e = 0-90 and P = 5 


le 


R-0 40 42 44 46 +8 06222 0-044 
38 +0 42 44 46 1-333 0556 0-167 
R-2 3-6 40 42 44 — 1556 0667 06370 0-222 | 0-133 
34 36 38 40 4:2 1778 O778 O444 0278 O-178 
R-5 30 3-2 34 36 38 R=5 2-222 1.000 «0389 0-267 
20 | 22 24 2-6 28 3333 «1-556 0-963 0-667 0-489 


Cake impurity F=02 R F-01\F=<02 F-03 
q 


be fore washi ng % % 


62-3 55-3 50-9 45-9 0 0 0 0 
R=1 61-6 56-3 51-9 46-9 R=1 0143 6-120 0-008 
R=2 61-5 60-7 57-3 52-0 47-9 R--3 251 0-332 O:298 0250 0-204 
R= 3 61-5 61-3 58-3 54-9 49-0 R-3 0-398 0-504 0-464 O-391 O321 
R=5 61-6 62-3 60-4 56-2 51-2 0691 O851 O714 0-582 
R 61-7 R 1-413 1-770 


Cake impurity F=02 F=03 Total yield F=061 F=02 


on it instead of the filter cake being submerged as Much more practical experience will have to be 
in the filtration proper. gained in a number of other filtration processes, 
before it can definitely be decided, whether or not 
any assumptions must be replaced by other, more 


Invariably values of «<4 were found for a drum 
submergence of about 180° and a wash section of 
about 130°. The values of x depend on the nature of convenient ones. 
the base material processed, by which the values of Acknowledgement is expressed to the management 
and are influenced. of N.V. De Bataajsche Petroleum Maatschappij for 

(b) The porosities of wax cakes are strongly permission to publish this paper. 
dependent on the nature of the base material and 
vary between 70% and 95%. NOTATION 

(c) The porosities of wax cakes after washing are A  ~ total surface area of a continuous filter 
normally a few per cent lower than before washing. 4.4 — surface area of a segment of a continuous filter 

Taking account of these facts, computed and F = solid content of base material (vol.) 
measured figures of cake impurity were in good © = fraction of original liquid obtained in primary fil- 
agreement with each other. Differences of not more trate (vol.) 
than 2-3°% were found. P = total volume of solvent used on base material 

Greater deviations from the theory can be expected RB volume of recirculation liquid on base material 
to occur when cracks are formed in the filter cake. &, ~ volume of original liquid in R | R 
R, 
8 


Then the basic wash equation will not hold any more volume of solvent in R — Re +h, 
and {(q,/q-) will be unfavourably affected. - volume of primary solvent used on base material 


| 
VOL. 
l 
1952 
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1—-F+S8 
F 
in slurry 


: ratio of volume of liquid to volume of solids 


filtration time of filter segment per cycle 

wash time of filter segment per cycle 

specific resistance of filter cake — resistance of filter 
cake per unit of surface area, when unit volume of 
filtrate has passed through it 

pressure difference at which filtration takes place 
volume of primary filtrate obtained in filtration 
time T 

volume of liquid in cake produced in time 7 — volume 
of pores in cake 

volume of wash liquid applied = volume of wash fil- 
trate obtained in wash time 7), 


f| qe | volume fraction of original liquid washed out from 
cake 

r resistance of filter cloth per unit of surface area 

T a+l 

ratio of cake resistance to cloth resistance at the 


x constant in basic wash equation = 
kq 
moment when washing is started 
pore content of cake = liquid content of cake (vol.) 
viscosity of primary filtrate at filtration tem- 
perature 
viscosity of wash solvent at filtration temperature 
mean viscosity of liquid in pores during washing 
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Summary—It is shown that the rheological behaviour of a fluid bed, which is supported by a gas dis- 
tributor, changes markedly along its height. These changes are connected in a qualitative way with the 
tendency of the gas to segregate from the dispersion of solid particles. 


Résumé—I1 est démontré que le comportement rhéologique d'un lit fluidifié, au-dessus d'un distribu- 
teur de gaz, change fort en fonction de sa hauteur. Ces changements dépendent qualitativement de la 
tendance du gaz & se ségréger hors de la dispersion de particules solides. 


Among others, Morse {1} has pointed out that a 
uniform dispersion of solid particles in a vertical gas 
stream is essentially unstable. A region of relatively 
low density tends to extract air from an adjoining 
region of high density. This segregation is counter- 
acted by the tendency of the solid particles to fill 
up the cavities. In large beds a balance between 
these two exists which results in gas pockets moving 
upwards by virtue of the solids falling down along 
their sides. The resulting high degree of agitation 
is one of the technical merits of industrial fluid bed 
reactors, as far as homogenization of temperature 
and transport of heat of reaction is concerned. 

In laboratory fluidization columns the bed is 
generally supported by a gauze or a piece of cloth. 
Immediately above this, the dispersion of the solid 
particles in the gas is fairly uniform. At greater 
heights segregation of gas proceeds and the system 
becomes more and more heterogeneous. At a certain 
height “‘slugs’’ may appear, i.e. gas bubbles of the 
same size as the tube diameter. MatTueson et al [2] 
in analogy to bubble formation in gas-liquid flow 


through narrow tubes, connected the tendency to- 
wards “‘slugging’’ with the “viscosity”’ of the fluid 
bed. A measure for this was the torque required to 
rotate a paddle of #4” ,” at 200r.p.m. in a 
internal diameter fluid bed column. 


As the transportation of heat in a fluid bed is 


mainly governed by the mobility of the solid particles 
and as this mobility will also manifest itself in the 
rheological behaviour, we tried to get some information 


on the latter property of such a system. 

Fig. | shows the experimental arrangement in 
which the rotational speed of a dumb-bell was meas- 
ured as a function of the torque applied to the rota- 
ting shaft. In order to minimize disturbances in the 
fluid bed and the interaction between the two spheres 
the rotational speed did not exceed 30r.p.m. Data 
were obtained at various superficial gas velocities (up), 
heights of measurement (h) and total expanded bed 
heights (/,,,) with the following materials: river 
sand 150-210 u, 105-150 uw, 75-105 u, a 50/50 % mixt- 
ure of coarse and fine sand, silica gel 105-150 yu, 
and Seotchlite glass beads 150-210 u and 300-420 y. 


Volt 
T 
k 
q 
q 
= 30 


The fluidizing agent was air. Pressure drop and 
pressure gradients were measured at the same time, 
but these results will not 
be discussed here. 


An arrangement like 
this does not permit quan- 
titative evaluation of the 
rheological behaviour such 
as a shearing stress rate 
of shear relationship. Thus 


the results can only be 
discussed in a qualitative 
way. As, moreover, the dif- 


ferent materials appeared 
to behave similarly, in the 
following we will restrict 


8-600m 


| fitter cloth 


an 
2 


ourselves to a few relevant 
facts. 


Air in 


Force per sphere (F) 


Fig. 1. Schematic diagram pyrene sphere velocity (v) 
and dimensions of experi- 

mental arrangement. Fig. 2 shows a typical set 
of data obtained at the 
bottom of a fluid bed at different values of the air 


velocity up, expressed in terms of uo,, the velocity 


es 
70 


River sand 105-150 
lexy =17-% om 
h=-3-90m 


5 10 cm/sec 75 


V-e 


Fig. 2. Force per sphere vs. linear sphere velocity at 
the bottom of a fluid bed. 


at which fluidization sets in. Apparently here the 
fluid bed behaves like a non-Newtonian liquid, giving 
a relatively high resistance to flow at small rates of 
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shear. These are the same features as observed with 
concentrated solutions of polymers and suspensions 
of solids. The curves can also be considered as inter- 
mediate between the behaviour of a fixed bed with 
a very high yield value and that of a dilute dispersion 
of solids in a fluid. In this respect the flattening out 
of the curves at higher rates is significant. 

For a further qualitative discussion we will use 
the force F at v=5 cm/sec, i.e. as a measure of 
viscosity. If Stokes’ law could be applied, F; = 5x 10° 
dynes would correspond to a viscosity of about 4000 ¢P. 


3 


ae 


- + 
River sand 
— 05-150 
---- 150-210 
Lexy 17—18-5 em 


7 2 


uot 

Fig. 3. Relation between “viscosity’’ and gas velocity 
at different heights in a fluid bed. 


Influence of air velocity (ug) on F, at 
different stirrer heights (h) 
Passing from a “quiet” to a “turbulent” state of 
the fluid bed, it is seen from Fig.3 that ¥;, first 
decreases rapidly to a low value and then tends to 
rise again. The behaviour in the former region may 
be compared with that of a liquid in which a relatively 
small increase of the average distance between the 
particles causes a marked decrease in viscosity. In 
a liquid the dilatation is determined by temperature, 
in a fluid bed by the gas velocity. At higher gas 
velocities, this trend is counteracted by the in- 
creasing turbulence or kinetic energy of the particles, 
which is brought about by the upward movement of 

gas pockets. 
The minimum value of F, is smallest at the bot- 
tom of the bed. This indicates that the degree of 
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turbulence increases in upward direction, i.e. that the 
gas bubbles grow on their way up. This was also ob- 
served experimentally by the fact that at a certain 
height the rotation of the dumb-bell became irregular 
because the size of the bubbles exceeded the size of the 
spheres. Beyond this height the interpretation of these 
measurements has more to do with statistics of the 
number and size of ascending gas bubbles than 
with flow properties of an apparently homogeneous 
substance. 


Dependence of F; on stirrer height (h) 
at high gas velocities 


that F; increases with h. 
with total bed heights of 


From Fig.3 it was seen 
However, in experiments 
30-40 em we found that 
bed F; decreased again 
accompanied with a more and more uneven rotation 
of the stirrer. This can be ascribed to the increasing 
size of the bubbles: once the stirrer has been struck 
by a gas pocket, it has more opportunity to gather 
speed, so that at a certain torque the average rota- 
tional speed becomes greater. Comparing the results 
of low and high beds under the same circumstances 
we had the impression that segregation of gas proceeds 
more rapidly in the latter cases, 

Although the above data are very complex and 
only a qualitative discussion could be given, we 


in the upper half of the 


considerably, which was 


thought it worth while to draw some attention to 
the fact that the resistance of a body which is dragged 
through a fluid bed is very sensitive to structural 
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changes in it. As the rheological behaviour varies, 
at least in vertical direction and probably also in a 
horizontal plane, the same might be the case with 
heat transfer coefficients in a small bed to a body 
suspended in it. Further it must be doubted on 
these grounds whether heat transfer measurements 
on a laboratory scale (column I.D. 5-30 cm) can be 
extrapolated for fluidization in technical dimensions, 
in which a different kind of solid-gas distribution exists. 


This conclusion can be stated in a more general 
way. Even under carefully specified conditions, no 
definite physical properties can be ascribed to a 
bed of fluidized solids because it is essentially non- 
homogeneous and non-uniform. In this respect it has 


nothing in common with the liquid state. 


The experimental part of this work has been 
performed by D. Meyer TiMMERMANN THyssen and 
J.J. Harr. 


NoraTION 
F = force per sphere (dynes) 
the same at v = 5 cm/sec 
height of measurement (cm) 
height of expanded bed (cm) 
superficial gas velocity (cm sec) 
the same at the beginning of fluidization (cm/sec) 
sphere velocity (cm/sec) 
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Summary— This study gives an investigation on the critical mass velocity at which a powdered solid starts to 
fluidize in an ascending flow of gas. From the measurements on closely sized semples and mixtures of 
different sizes of carborundum, iron oxide (Fe,O,) and coke, using air, argon, carbon dioxide, nitrogen- 
hydrogen mixtures, town gas and methane as fluidizing gas, the proportionality of the critical mass velo- 
eity Gy, with the reciprocal value of the kinematic viscosity 4/@ of the gas and with the bed density at 
maximum porosity @),, and the square of the particle diameter d was confirmed in the region Re < 5. 
By the introduction of the so called effective diameter d, (which is the diameter of the sphere having the 
same number of particles per unit volume of bed at maximum porosity) and of a generalized shape factor B, 
which includes the usual porosity term (‘‘Kozeny factor’’) it was demonstrated that instead of the porosity 
and particle diameter, the principal property which determines the resistance to gas flow is the number of 
particles per unit volume of bed. These considerations resulted in an equation for the critical mass velocity 


reading 
Re, = 


0,00123 2 
B 
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in which when substituting for d the effective diameter d, the generalized shape factor B, equals, even for 
definitely non-spherical particles like carborundum and ironoxide (see Fig. 2a, b and c), practically unity. 

After the introduction of the packing ratio m, being the ratio of bed density and bed density at maximum 
porosity a modified “Carman-Kozeny” equation for the pressure loss in a fixed bed is suggested where the 
usual porosity term is substituted by a power of the packing ratio. 


Résumé—L étude porte sur un examen de la vitesse critique A laquelle une substance solide en poudre 
commence & se fluidifier dans un courant de gaz montant. Les mesurages effectués sur des échantillons 
de grosseur exacte ainsi que sur des mélanges de carborundum, d’oxyde de fer (Fe,O,) et de coke, différant 
en grosseur, alors que lair, argon, le bioxyde de carbone, les mélanges azote-hydrogéne, le gaz de ville et 
le méthane servent de gaz fluidifiant, ont prouvé le rapport entre la vitesse critique de la masse Gy, les 
valeurs réciproques de la viscosité cinématique “0 du gaz, la densité du lit & porosité maximum 9, ,, et le 
carré du diamétre de particules dans la région Re <5. 

En introduisant le diamétre effectif d, (c'est & dire le diamétre de la sphére contenant la méme quantité 
de particules par unité de volume du lit & porosité maximum) et le facteur de forme genéralisée B renfermant 
le terme usuel pour la porgsité («facteur Kozeny») il a été démontré qu'au lieu de la porosité et du diamétre 
des particules, la principale qualité déterminant la résistance par rapport au courant de gaz, est la quan- 
tité de particules par unité de volume du lit. De ces observations, il résulte équation pour la vitesse cri- 
tique de la masse: 

0.00123 0%, 
B 


dans laquelle, si d est substitué par le diamétre effectif d,, le facteur de forme généralisée B, est égal a 1, 
méme en cas d'emploi de particules qui n'ont pas la forme sphérique tels que le carborundum et oxyde de 
fer (voir figures 2a, b et c). 

Aprés lintroduction de l'indice de la “packing ratio’ m, c'est a dire le rapport entre la densité du lit 
et la densité du lit & porosité maximum, il a été proposé un rapport Carman-Kozeny modifié pour la perte 
de pression d'un lit fixe dans lequel le terme usuel pour la porosité a été substitué par une puissance de 
Vindice de la “packing ratio.” 


INTRODUCTION extensive study on some features of the fluidized 


During the last ten years the technique of fluidization ‘state. The first object, which is dealt with in this 


has found a rapidly growing application in the chemi- 
eal industry. The main characteristics of this unit- 
operation which make its application so attractive, are: 

1. the ease, by which large masses of solids can 
be transported from one reactor to another, per- 
mitting in many cases an efficient utilization of heat 
and a ready development of a continuous process; 

2. the rapid mixing of the solid particles, resulting, 
even in the case of strongly exothermic or endothermic 
reactions, in a very homogeneous temperature 
throughout the reactor and in high rates of heat 
transfer to heating or cooling surfaces; 

3. the possibility of using finely divided solids 
without suffering an abnormally high pressure loss 
at reasonably high space velocities. 

The main impetus to the development of the 
fluidization technique was given by the oil cracking 
industry, where an accurate control of the reaction 
temperature and an easy regeneration of the solid 
catalyst combined with an efficient utilization of 
the heat of regeneration is most important. 

Although much work has already been done in 
the field of fluidization, there is still a great lack of 
fundamental knowledge as regards the physical 
properties of the fluidized state of gas-solid systems. 
This consideration induced the authors to start an 


paper, was to study the dependency of the critical 
mass velocity at which fluidization starts on the 
physical properties of gas and solid. The second 
object—-which will be dealt with in a following 
paper—was a study of heat transfer characteristics 


under widely varying conditions of gas and solid. 


THEORETICAL CONSIDERATIONS 


When a gas flows in upward direction through a bed 
of powdered solids, it suffers a pressure loss, which as 
a function of the mass velocity is schematically 
shown in Figure 1. 

For low flow rates the pressure loss increases 
linearly with the mass velocity until a point is reached 
at which the pressure loss equals the static weight 
per unit cross sectional area of the solid powder, 
Ai a further increase of the rate of flow the bed of 
solids expands until the particles are free to move 
with respect to each other. From this moment the 
pressure loss is substantially constant and equals, 
independently of the mass velocity the static pressure 
of the solid. This can be expressed by the following 
formula * 


(1) 


* See nomenclature at the end of the paper. 


Volt 
No, 1 — 1951 
or, when 0 is neglected against 9,, by: 


Ap=go,l—o)L=go,L. (2) 


This relation has been confirmed by many inves- 
tigators. Deviations are mostly due to non-essential 
complications like channelling, electrostatic effects or 
frictional forces of the wall, when using smaller tubes. 

It seems plausible to define the critical mass 
velocity Gg at which fluidization starts, as the mass 
velocity belonging to the point of intersection A of 
the two straight lines in Figure 1. However for a 
given powder this point ist not unambiguously deter- 
mined, as the slope of the line for the fixed bed region 
depends on the original packing of the powder. 

This difficulty is avoided by measuring the 
pressure loss-mass velocity relation with decreasing 
flow rate, starting in the fluidization region. In this 
case the fixed bed is settled in a reproducible way at, 
what we will call maximum porosity. 

This means that we define the critical mass velo- 
city G, for fluidization, as that mass velocity, at 
which the pressure gradient in the fixed bed under 
the condition of maximum porosity equals the static 
weight gradient of the solid. 

The first authors, who gave this definition of the 
starting point of fluidization, were Leva, Grummer, 
Werntraves and Poticurk [7] who introduced the 
concept of “minimum fluid voidage” which is vir- 
tually the same as the maximum porosity of the 
fixed bed. 

The starting point of fluidization can also be 
measured at an increasing flow rate when owing to a 
preceding fluidization, followed by slowly reducing 
the gas velocity down to zero the packing is of maxi- 
mum porosity. In this case the above mentioned 
expanding of the bed, when passing at increasing 
flow rate the critical mass velocity, practically does 
not occur, 

In order to calculate the critical mass velocity Gy 
for a given gas-solid combination the relation between 
pressure loss and gas flow rate in a fixed bed must 
be known. This is a very old problem, which was 
firstly discussed in a publication by p’Arcy [1]. Espe- 
cially the work of Biake [2], Kozeny [6] and Car- 
MAN [3] should be mentioned. BLake introduced the 
concept of specific surface per unit volume of packed 
bed a, which led to the following equation: 

k 


_ Ap me 
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This relation, which is only valid in the region of 
viscous flow, i.e. in the region where the pressure 
drop is proportional to the mass velocity, is con- 
firmed by experiments, as was shown by CarMAN [3]. 
who collected the available data from literature. For 
definite geometric particle shapes, from which the 
surface area could be calculated Carman always 
found approximately the same constant k~5. He 
therefore postulated that k= 5-0 for all shapes and 
considered relation (2) as a possibility for determining 
the specific surface of powders from pressure drop 


measurements. From later experiments, however, 


(e.g. OMAN and Watson [10}) it appeared that when 
3 A 


Fluidization region 


| 


& Mass velocity 
Fig. 1. Pressure drop—flow rate relation. 
the particles stronsiy deviate from the spherical 
shape, k can attain much higher values (Raschig 
rings k 20). 
6(1 
eq.(3) is transformed 


For spheres, where a = , : 


into: 


Ap od Cc 


The porosity term is usually called the 


(1 — 
Kozeny-factor. 

From the Carman value k=5 follows C=—90 
while Leva, GrumMeER, and PoLicnik [7] 
deduce from their experiments C = 100. We prefer the 
result C =77 which was obtained in the excellent 
experiments which Lewis, GrLLILaAND and Baver [8] 
performed with glass spheres of very homogeneous 
diameters. 

After the diameter of non spherical particles has 
somehow been defined, eq. (4) is formally written as: 

_Ap ed 


f=) (5) 


The shape factor 8 (> 1) is a measure of the deviation 
of the particles from the spherical shape (f= 1). 
The shape factor introduced in this way must be 
considered as an empirical constant, which can only 
be determined from a measurement of the pressure 
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loss, its absolute value depending on the more or 
less arbitrarily chosen diameter such as nominal and 
statistical diameter or the diameter deduced from 
elutriation-, sedimentation or sieve analysis (see 
DALAVALLE [4}). 

From eq. (2) and (5) the critical mass velocity can 
be calculated by eliminating 4p and putting = 77 
and 6=46,. This results in the following relation, 
firstly derived by Leva c.s. {7}. 
q 1 


(6) 


where g = 


_ a dimensionless number introduced 


among others by Witnetm and Kwavuk [12}. 

It was the object of our experiments performed on 
different solids of varying particle size, in conjunction 
with different gases, to test the usefulness of eq. (6) 
for predicting the critical mass velocity for a given 
powder. Because eq. (6) is only valid in the region 
of viscous flow, we restricted our experiments to such 
conditions, that the pressure drop in the fixed bed was 
always proportional to the gas flow rate. The limit of 
the viscous region is given by the condition Re < ~ 10, 


EXPERIMENTAL 
Apparatus 


The pressure loss-flow rate relation according to 
Figure | was measured in a very simple apparatus, 
consisting of a glass tube of 8-5 cm in diameter having 
at its lower end a removable conical gas inlet made 
of brass which was provided with a sintered glass 
filter. A box made of thin brass plate was placed on 
the upper end of the glass cylinder which for reason 
of its wide diameter (17 cm) was used as dust trap. 
The manometer connections were mounted on the 
conical inlet and the dust collector, the pressure drop 
was measured with a water manometer and corrected 
for the pressure loss in the empty tube. 

The height of the settled bed was always about 
20 cm. Some experiments made with different heights 
of bed proved that the resulting value of Gy was 
independent of the height of bed. 

The bed density at maximum porosity was deter- 
mined in a separate all-glass apparatus, consisting of 
a cylinder ( © 6,2cm) provided with a conical inlet and 
a sintered glass filter. A known weight of powder was 
brought into the cylinder, fluidized, whereafter the 
flow rate was slowly reduced to zero. The volume of 
the powder was then read from a calibrated scale on 
the cylinder wall. 
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Data of used materials 
The solid materials used in the experiments are 
mentioned in Tables and 2. Table l shows the 
homogeneous fractions, Table 2 mixtures of different 
sizes, which were prepared by mixing known weights 
of homogeneous fractions (the sieve analysis is given 


* in Table 3). 


Table 1. Properties of homogeneous fractions 


Sieve- 
limits (u) 


75% 82 97 3180) 1680 
90100 113 
110-125 
150-175 
175-210 
210-240 
75-110 
75-110 
125-175 
150-210 
210-240 
240-300 
300-420 
420-500 
500-800 


carborundum 


iron oxide 


(Fe,O,) 75-110 


Table 2. Mixtures 


Sreve-limits 


(a) d(u) d,(u) 


carborundum 
carberundum 


75-210 128) 1663) «(0-477 
50% -75-90 
50%-210-240 120 1825) 0-426 
carborundum 50-210 1700) 0-465 
coke 5-210 70) 0-564 
iron oxide 


(Fe,0,) 50-300 2870 


As diameter of the homogeneous fractions is chosen 
the arithmetical mean of the sieve limits; the diameter 
of mixtures a and + is calculated from the weight 
fractions x, of the components of diameter d, from the 


relation: 


1 
dwixt xi/d; (7) 


The reason why we made this choice out of a 
number of other possibilities will be explained later. 


VOL. 
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0-471 
0-496 
O-487 
476 
0-475 
coke 0-564 
0-595 
0-600 
0-624 
0-631 
5180 2620 0-494 
sud - 
a 
b 
d 
€ 
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Table 3. Sieve analysis of mixtures in °o 


Sieve 50- 65- 
on 
65 75 


100 


Sclid 


carborundum 20-1 13-8 
carborundum 
 earborundum 9-6 75 2:7 23-3 
d coke 2-3 3-2 11-2 14-1 
iron oxide 
(Fe,0,) 44 46 12-4 12-1 


100- 125- 150- 175- 210- 240- 275- 
125 150 175 210 240 275 300 


19-7 12-5 12-5 12-4 
15-0 10-0 10-0 
24-5 12-2 1-0 1s 
50 64 12-4 4-2 9-9 


The mixtures c, d and ¢ have a lower sieve limit 
of 50u. Since a reliable sieve analysis is difficult to 
carry out in this region, a sieve diameter of these 
mixtures is not mentioned*. The meaning of the 
effective diameter will be explained later. 

The bed density at maximum porosity o,,, was 
determined in the aforementioned apparatus, the 
material density 0, was determined pyknometrically. 

The coke particles, when viewed under a micro- 
scope show relatively big holes and cavities in the 
surface (see Figure 2e). When determining the mer- 
cury-density as well as the helium-density of two 
extreme sizes, we found for the 420-500u fraction 
Oy, 1990 kg 1800 kg'm® and for the 
fraction oy, 2040 kg a4, — 1800 kg 
Apparently both o,, and oy, are practically inde- 
pendent of particle size in this region. 

For further calculations the mercury density was 
used as it is quite clear that the pores accessible for 
helium but not accessible for mercury (pore diameter 
«10y) hardly contribute to the passage of gas flow 
through the bed. 

From the relation =o, (1 


ean be calculated. 


the maximum 
It is clear, that the 
high porosity of the coke beds, whith increases still 
further at increasing particle size, is due to the porous 


porosity 


m 


structure of the particles. It is probable that a large 
part of this porosity will not contribute to the passage 
of the gas flow. For these kinds of materials it is 
impossible to distinguish between particle- and bed 
porosity. 

In contradistinction to the results of Leva c.s. 
who found that the porosity of sand increases with 
particle size, the porosities of our carborundum 
samples show no systematical relation of this kind. 

Figure 2 gives the photomicrographs of some of the 
materials used. 

* Since in eq. (7) the weighted mean of the reciprocal 


diameters is calculated, the influence of the particles with 
smallest size on the ultimate value of dpjx¢ is highest. 
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The gases employed in our experiments are men- 
tioned in Table 4. The density is taken at 25°C and 
Latm, which are the normal experimental conditions. 


Corrections for temperature and barometric deviations 


Table 4. Properties of gases 


Gas (25°) 10° 0 wo. 10° 
98% H,+2% N, 0-93 0-104 8-95 
80% H, + 20% N, 1-31 0-297 4-41 
65% H, + 35% N, 1-495 0-458 32 
45% H, + 55% N, 1-64 0-671 2-44 
25% H, + 75% N, 171 0-886 1-93 
Air 1-85 1-185 1-56 
Argon 2-23 1-65 1-35 
Carbon dioxide 1-50 1-82 0-825 
Towngas 1-22 0-349 3-50 
Methane (raw) 1-24 0-700 1-77 


and for the pressure loss in the bed were not made 
since the collective deviations were always less than 
2%. The viscosity was measured against air as a 
standard in an apparatus according to RANKINE {11}. 
A comparison between our results for CO, and Ar and 
the values taken from literature proves the reliability 
of the method. 

The measured values of town gas and raw methane 
were in good agreement with Manne formula for the 
viscosity of gas mixtures, reading: 


% 


8 
*mixt 
The results found for the N,-H, mixtures proved 

to deviate up to 9% from eq. (8). 


Experimental results 
The results of the Gy determination for the different 
materials and gases have been compiled in Tables 5 
and 6, the measured pressure loss-flow rate relations 
are shown in Figure 3a-f. 
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Fig. 2e. Coke 500-800. Magn. 30>. 


Fig. 2a-e. Photomicrographs of the materials use«| 
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Table 5. Results on homogeneous fractions 


Solid Was G Re B B LB 
carborundum 7500 air O-0504 31-6 1-30 0-66 0-92 
80-100 air 0-376 O-OL7S 49-0 1-31 0-67 O-94 
110-125 air (484 0-0260 0-165 1-68 0-66 O-83 
air 441 0-054 472 236 1-45 0-62 
175-210 air 0-392 0-069 716 402 1-43 0-69 0-99 
210 240 air O-OS4 1-023 647 1-60 78 1-03 
75-110 air (376 OOLTS (0880 47-6 
75 110 H, 0-376 0-00266 00275 15-8 
75-110 H, 80-00% 0-376 0-0059 0-0423 23-6 
75-110 H, 61-09% O-376 0-0079 0-0488 29-0 — 
75-110 H, 45-6% O-376 -O103 O-0591 33-9 
1-37 0-69 
75-110 H, 23-5% 0-376 O-O135 0-0737 41-9 
75-110 Ar 0-376 0-0202 O-O851 45-5 
75-110 CO, 0-376 0-032 0-201 110.0 
75-110 towngas (376 0-0071 00546 32-1 
75-110 methane O-376 O-O155 62-2 
coke 75-110 air O-919 0-0146 0-0733 21-6 
75-110 H, 0-0232 7-2 
75-110 H, 80-8 O-0045 0-0322 mean value 
75-110 H, O-0065 0-0403 12-9 
1-90 39 0-60 
75-110 H, 23-9 O-919 O-OLLS O-0615 19-3 
75-110 Ar (0166 0.0693 20-7 
75-110 CO, 0-0237 0-147 
125-175 air 1-174 0-0295 0-237 83-5 2-68 O43 0-64 
150-210 air 1-066 O-O41 0-401 153 2-59 0-47 O55 
210-240 air 1-286 0-067 O-S15 287 2-83 0-43 0-56 
240-300 air 1-346 0-100 1-488 509 2-99 0-42 0-65 
300-420 air 1-718 0-161 3-14 1070 B81 0-42 0-56 
4206-500 air 1-85 225 5-65 2270 4-52 0-65 
500-800 air 240 0-350 12-3 5810 7°36 0-58 0-60 
iron oxide 75-110 air O47 0-028S8 0-145 69-5 
75-10 H, O47 O-0048 24-2 
75-110 H, 0-47 0-0097 0-0686 33-6 mean value 
75-10 H, O47 O-0858 41-6 147 O59 
75-110 H, 47 O-O179 L004 40-8 
75-110 H, 241% O47 00-0249 61-6 
75-110 Ar O47 O-0356 O-147 66-6 
75-110 CO, O47 O-0545 O-334 162 
Table 6. Results on mixtures 
lieve 
Solid Ge Gy Re, 4 B B B, 
a carborundum 75-210 air 0-396 0-025 0-149 75-2 1-30 
carborundum 50-75-00 air 0-235 O-O194 0-127 107-2 1-29 1-04 
210-240 air 
carborundum 5-210 air O35 0-0193 O-118* 83-6* — — 
d coke 5-210 air 1-346 O-O195 0-136* §2-2* 0-47 
iron oxide 50-300 air 0-29 0-033  0-248* 214* 1-06 


* Calculated with the effective diameter d,. 
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Corborundum- gases 
bad 98,5% Hy + 1,5 %N, 
© +200 % Ny | 
61,9 + 38,7 %N, | 
46,8 Vo, + 548 %N, 
23,5 + 76S | 
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56789107 
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Fig. 3a. Pressure drop—mass velocity relation for carbo- 


rundum 75-110u with different gases. 
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Fig. 3c. Pressure drop—mass velocity relation for mixtures 


of Table 6. 
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Fig. 3e. Pressure drop—mass velocity relation for iron 


oxide 75-1104 with different gases. 
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Fig. 3b. Pressure drop— mass velocity relation for carbo- 

rundum fractions of different particle size. 
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Discussion oF THE RESULTS This illustrates the dependency of the shape factor 
Application of the Leva-Carman eq. (6) on the method in which the bed porosity is determined. 
In Figure 4 the experimental values of Gy for the same 
fractions of carborundum, iron oxide (FegQ,4) and 6 0 oy 20% 
of the different gases. The points are correlated with Si.’ pS 
a reasonable accuracy by a straight line with slope -1, . -» Carbon dioxide 
proving that eq. (6) accounts in the right way for the 
© Methane (raw) 


properties of the gas. 


In Figure 5a 


carborundum data with air. The values of the shape 


plot is given for all 


+ 


factor f for fractions of different particle sizes, mentio- 
ned in Table 5, show a maximum deviation of 15% 


from their mean value 1-45 proving that the shape 
factor of carborundum is reasonably independent of : 
the particle size. 
Figure 5 also gives the result found for the carbo- 4 
rundum mixtures a and (Table 6). Both mixtures msec 
correspond to four points, lying on a straight line Fig. 4. The critical mass velocity Gyvs the kinematic 


. viscosity “0. 
with slope +13 and being calculated from four dees 


different definitions of the mean diameter (see legend 


Moreover, the problem with which we are con- 
Figure 5). Other possible definitions of the mean 


fronted is probably of a more general character and 
diameter of a mixture are outside the range of the 


four chosen (see DALLAVALLE [4}). eo 


| 
+ 


From the position of the different points, espe- 
cially those of mixture 4, it is quite clear, that only the 


— 


definition correlates the results of mixtures with those 


of homogeneous fractions. | 


‘ 
The caleulated values of the different coke i: 

fractions prove to increase with particle size. The ae a a5 
reason is obvious: according as the particle size in- ‘ ° 150-176 ja 

creases an increasing part of the determined porosity 
3 
is to be considered as particle porosity this Mit 
j , é 150% —75— 90 
contributes only to a small extent to the passage of be 
the gas flow. The consequence is, that the resistance ) Eard 

of the bed is much higher than should be expected 
from the high porosity values, resulting in a high # 97 fs 


value which increases with particle size. 53 
Fig. 5. Rey plot of the carborundum fractions. 
m 


It might be reasoned that # is only an empirical 
constant, so that abnormally high values of more or Illustration of the correctness of the d= v 
less porous particles do not violate the correctness of a ld 
eq. (6). However if we had used the helium—instead 
of the mercury-density for calculating the porosity of | not only restricted to porous particles in the limited 
the coke beds the resulting shape factor would have — sense of the word. The aforementioned high values 
been still higher. of k in the general eq. (3) for particles such as Raschig 


definition for mixtures. 
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rings show the same tendency. High porosities, when 
caused by peculiar particle shape, tend to give high 
shape factors. 

As was pointed out already by Happens [5] it is 
very probable that the porosity term 


in eq. (5) only accounts for a change in porosity caused 
by varying the density of packing, while the influence 
of an increase in bed porosity by altering the shape 
of the particles is much less. 


Generalized jorm of eq.(6) 
As was argued in the preceding section it is probably 
wrong to account for the influence of bed porosity 
on the resistance of gas flow in one definite porosity 
factor for all kinds of particles. 

Especially for the application of the resistance 
law in packed beds to the starting point of fluidization 
formulated in eq. (6) a generalization can be made 
which eliminates the porosity term. 

The maximum porosity of a fixed bed, +. the 
porosity of the loosest stable packing of particles 
resting by contact on each other, is completely 
determined by the geometrical shape of the particles. 
This is confirmed by some experimental results 
obtained with spheres. In our experiment performed 
with glass spheres of | mm diameter we found a 
maximum porosity of 0-406. Lewis, GILLILAND and 
Baver [8| mention as the highest porosity they could 
attain for glass beads of 150u a value of 0-409, 
Mickiey and Tritiine [9| determined the bulk 
density of beds of Scotchlite glass beads of homo- 
geneous diameter varying from 70 to 4504. From 
these data the porosity, which in this case must 
be expected to be somewhat smaller than the maxi- 
mum porosity, could be calculated, the resulting 
values varying only from 0-390 to 0-397. 

These results indicate that the maximum porosity 
d of a fixed bed of spheres of homogeneous 


diameter has a closely determined value which is 


m sph 


independent of particle size. Hereafter we postulate 
sph 0406. 

Now if we assume that for all particle shapes 
+, is completely determined by the geometric proper- 
ties of the particles, it is possible to combine the pro- 
duet of the porosity term and the shape factor f in 
eq. (6) into one generalized shape factor. Moreover, if 
a quantitative value of the bed porosity makes no 
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sense, it is also necessary to reject the concept of 


particle density and to use only the well defined bed 
density at maximum porosity. 
In this way eq. (6) can be generalized to: 


Reg Xq (9) 


ge 


where is the same number as in 
eq. (6) after substituting o, by o,.,,. 
For spheres X can be interpreted from eq. (6) as: 


Xn = 23. (1) 
(i sph)* 
Eq.(9) can now be written as: 
Rey 040123 13 at) 


B 


where the generalized shape factor B equals unity 
for spheres. 

The calculated B-values mentioned in Tables 5 
> 1, Bon the 
This means, that when the 


and 6 show that, whereas / is always 
contrary is always < |. 
resistance to gas flow of a bed of particles of various 
shapes is compared with the resistance of a bed of 
spheres of the same diameter, the resistance of the 
former will be higher than that of the latter when 
the comparison is made at equal porosity, but lower 
when said resistances are compared at maximum 
porosity. 

The spread and the deviation from unity of the 
B-values for the different carborundum fractions 
(mean value 0-68, mean deviation 5°) is practically 
the same as for the p-values (mean value 1-46, mean 
deviation 8°). This proves that, apart from the 
considerable simplification obtained by the elimination 
of the porosity term, the introduction of the generalized 
shape factor for these nonporous particles is at least 
equally satisfactory. 

For the different coke fractions the improvement 
is remarkable, the variation in B (0-39-0-58) being 
very much smaller than in # (1-9-7-4). 

For mixtures the elimination of the bed porosity 
seems to be less suitable as is demonstrated by the 
result of the carborundum mixtures a and + in Table 6. 
Mixture 6 has a much lower porosity 4,,, showing 
simultaneously a higher B-value. The carborundum 
mixture a, however, having a maximum porosity of 
the same order as the homogeneous fractions, shows 
the same value of B. Moreover the f-value for both 
mixtures is the same as for the homogeneous fractions. 
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On account hereof it may be suggested that after 
elimination of the porosity, the influence of a decrease 
in porosity, effected by mixing components of dif- 
ferent particle size, is inadequately accounted for. 

In the following section it will be shown, however, 
that it is possible to rule out this discrepancy by 
introducing the proper diameter definition, still without 
making use of the numerical concept of bed porosity. 


Effective diameter 


As was pointed out before it is rather unsatisfactory, 
that the shape factor (both B and #) depends on the 
arbitrarily chosen definition of the diameter of the 
particles (in our case the sieve diameter). There is 
however a diameter definition which is completely 
adapted to the present problem. 

The definition of this diameter, hereafter called 
the effective diameter, reads as follows: the effective 
particle diameter d, of a powdered solid is the diameter 
of the sphere with the same number of particles per unit 
volume of packed bed, both counted at maximum poro- 
sity. 

This definition is analogous to the so-called 
“nominal diameter,” the only difference being that 
whereas the “‘nominal diameter” is based on the 
number of particles per unit volume of solid, the 
effective diameter is based on the number of particles 
per unit volume of bed at maximum porosity. 

For mixtures the definition of the effective dia- 
meter is also directly applicable. 

The effective diameter d of all materials was 
determined by counting the number of particles in 
a sample of known weight. The number x,, of particles 
per unit volume of bed at maximum porosity then 
follows directly from the bed density 9,,,,. 

From the maximum porosity 4,, .,,, = 0-406 of 
a bed of spheres it follows that the number x,, .., 


of spheres of diameter d per unit volume of bed at 


maximum porosity amounts to 0-594 —s . Putting 


n, equal to gives: 
— 1088 


Determining the number of particles in a known 
weight of powder could always be done by spreading 
the sample on a piece of graphical paper and counting 
every square centimeter separately. The weight of 
the samples varied from 0-0036 g to 0-17 g, the highest 
number counted amounted to 7200. 


(12) 
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All samples were taken by reduction of the whole 
quantity of the material necessary for the pressure 
drop measurement to a countable number of particles 
with the aid of a small type riffle sampler. 


The resulting values B, (Tables 5 and 6) of the 
generalized shape factor calculated on the basis of 
the effective diameter are remarkable. The values 
of the homogeneous carborundum and iron oxide 
fractions are very near to unity, showing that with 
the concept of a generalized shape factor and the 
introduction of the effective diameter, these particles 
behave almost like spheres. 


This means, that the principal magnitude which 
determines the resistance of a bed at maximum poro- 
sity is, irrespective of the numerical value of the 
porosity, the number of particles per unit volume 
of bed. 

Even the carborundum mixture 6 with its rela- 
tively low porosity value now shows a shape factor 
B,=1-0, which deviates only 7% from the mean 
value 0-94 of the homogeneous fractions. 

The B,-values of the different coke fractions show 
a mean deviation of only 6% from their mean value 
0-60 and are also much nearer to unity than the B- 
and #-values. 

Besides the two carborundum mixtures mentioned 
before Table 6 also shows the results found for three 
mixtures with a lower sieve limit of 50u. Since a 
reliable sieve analysis is difficult to carry out in this 
region only the effective diameter was determined, 
the resulting values of B, being also quite satis- 
factory. 

The elimination of the porosity term and the 
introduction of the generalized shape factor and 
effective diameter definitely improves the correlation 
of particles of different shape and structure. Espe- 
cially for “porous”’ particles like coke, for which a 
numerical concept of bed porosity does not make 
any sense, it leads to acceptable values of the shape 
factor. 

Moreover substantially non spherical particles like 
carborundum and iron oxide (see Fig. 2) which must 
be considered as normal representatives of a large 
class of powdered solids obtained by crushing and 
grinding, show a generalized shape factor ~1. 


This means, that after determination of the effee- 
tive diameter the starting point of fluidization can 
in many cases be predicted with reasonable accuracy 
from eq. (11) by putting B,=1. 
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The bulk density, normally determined by filling 
a graduated glass with a loosely packed material, 
deviates only to a small extent from the bed density 
at maximum porosity, #.¢. the porosity after fluidiza- 
tion and subsequent slow reduction of the flow rate 
down to zero. Because according to eq. (12) the effec. 
tive diameter is only proportional to the cube root 


Corborundum 90-100 
150-175 
210-200 
75-110 4 
270-200 
500-800 
ronoxide 75-110 
Glass spheres 150 wiLewss cs) 
7 
105 495 120 425 130 
Packing ratio m 


Fig. 6a. Ratio of resistance at minimum and maximum 
porosity vs the packing ratio m. 


Resistance ratio 


|)” Legend see fig. €a 


; ++ +4 


7? 2 3 
9-6 


Fig. 6b. Resistance ratio es the ratio of the 
Kozenyfactors. 


of the packing density, the normally determined 
bulk density is quite sufficient for calculating the 
effective diameter. 


The influence of the packing density on pressure drop 


The conclusions stated in the preceding sections were 
all based on the measurement of pressure drop in a 
fixed bed at maximum porosity. It remains necessary, 
however, to formulate the influence on pressure drop 
of a change in porosity by variation of the packing 
density. It is a logical consequence of our preceding 
observations to formulate this influence also on the 
basis of the number of particles per unit volume 
of bed. 


At maximum porosity the pressure drop-flow rate 
relation can be derived from eq.(11) and (2) giving: 


j Ap ode 1 410 


L 2@ Re’ (13) 


In order to account for the increase of pressure 
drop when making a denser packing, we introduce 
the packing ratio m, which is the ratio of the number 
of particles per unit volume of bed at the considered 
porosity and at maximum porosity, which is the same 
as the ratio of bed densityand bed density at maximum 
porosity. 

This can be expressed by the formula: 


m (14) 


Ob m 


An extension of eq.(13) to other porosities can 
now be written as: 


Ip od, 410 


L p> (15) 


where y(m) is a function of the packing ratio m, 
with w(l)=1. 

In order to investigate this assumption the pressure 
drop-flow rate relation was measured for some ma- 
terials at the highest density, that could be obtained 
by slowly filling the apparatus, tapping simultaneously 
on the tube wall (see Table 7). 


Table 7. Dense packings 


Sieve 
Selid limits (u) * 


carborundum 100 1855) 214 1-10 
earborundum 150-175 O208 114 
exrborundum 210-240 1845 0420 0-220 1-10 
coke 75-110 922. 0-487 0439 (1-16 
coke 500-800 7 0-610 1495 
iron oxide 75-110 2845 O-441 1-10 


Intermediate values of the porosity were not in- 
vestigated, because it must be seriously doubted, 
whether it is possible in this case to obtain a homo- 
geneous porosity throughout the whole bed. 

The results of these experiments are given graphi- 
cally in Figure 6a and b. In these figures, the ratio 
of the pressure drop at the densest packing and the 
pressure drop at maximum porosity both at the same 
height of bed and mass velocity is given as a function 
of the packing ratio resp. of the ratio of the Kozeny 
factors for the two porosities. 


VOL, 


1952 


+ +—+—+ r+ 
48 


Vol.T 
No. 1 — 1951 


These figures also show the ratio of the results 
of Lewis, GrLuitanp and Baver [8] for the two 
extreme porosities of a bed of glass spheres of 150 u 
diameter. 

In Figure 6b the point of Lewis c.s. for spheres 
lies, as should be expected exactly on the straight line 
according to eq. (4). Our own results obtained with 
non spherical particles, however, deviate considerably 
from this line. In Figure 6 a all points are reasonably 
well correlated by a straight line with slope + 6-8. 

Here too it appears that the pressure drop in a 
packed bed of particles of various shape is better 
correlated by the introduction of the number of par- 
ticles per unit volume of bed than by the use of a 
porosity term. 

As a final equation we suggest instead of eq. (5): 
L 2G? \m/! B, 


410 


i Re,’ 


(16) 


A further extension of the experimental data is 
necessary to test this formula on its ultimate value. 


Remark on the maximum porosity of a bed of spheres 
There are four systematic packings of spheres [4] 
0-260), 
0-476) and two 
intermediate arrangements, the tetragonal spheroidal 
(d= 0-302) and the orthorombic packing (0 = 0-395). 
It is interesting to remark, that the experimental 


the closest ‘‘ rhombohedral” type of packing (6 
the most open cubical packing (6 


value of the maximum porosity of a bed of spheres 
-0-406 deviates very little from the calculated 
porosity of the orthorombic packing. The agreement 


maph 


with the bulk porosities determined by MickLey and 
[9}, which varied from 0-390-0-397, is even 
still better. This agreement seems to be hardly acci- 
dental, and it may be supposed, that the orthorombic 
packing must be considered as the loosest stable 
arrangement. 

In this connection it might be of interest to study 
by a geometric statistical analysis the random packing 
of spheres, obtained by arbitrarily adding the particles 
one by one. 


NOTATION 
a — specific surface per unit volume of packed bed m*/m® 
(— 0-305 ft2/ft®) 
B = generalized shape factor calculated with sieve dia- 
meter 
B, = generalized shape factor calculated with effective 
diameter 
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d = sieve diameter m (-- 3-28 ft) 
d, -- effective diameter 


d, — diameter of a component in mixture 
modified friction factor 
acceleration of gravity 9-8 m/sec® (— 4-13 
10° ft/hr?) 
G = mass velocity in empty tube kg/m? sec 
( 738 Ib/ft® h) 
Gy — critical mass velocity where fluidization starts 
L — height of unexpanded bed m( = 3-28 ft) 
m packing ratio 
n- number of particles per unit volume of bed 
n,, tumber of particles per unit volume of bed at 
max, porosity 


msph ~ “umber of spheres of equal diameter per unit 
vol. of bed at max. porosity. 
Gd 

Re = Reynolds number 

Gd 

Reg 
Gd, 

Re, 


x; = weight fraction of a component in mixture 
shape factor according to the Leva-Carman eq. 
bed porosity 


maximum bed porosity 
Om 0406 maximum porosity of a bed of spheres 
of equal diameter 
Ap pressure drop V/m*( 8-7 10° lb force/ft*) 
— dynamic viscosity kg/m sec( 24201b/hr ft) 
vy kinematic viscosity m*/sec (= 3-87 10*ft®/hr) 
density of gas kg m* ( — 0-0624 Ib/ft*) 
v, ~ density of solid 
density of bed 
0, — ved density at maximum porosity 
Clb m 
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R. C. L. Bosworrn; Physies in Chemical Industry. MacMillan 
and Co., Ltd., London, 1950. XX + 928 pp, 15 « 22 em, Price 
bound £ 3-10s 


In his introduction the author says that by writing this 
book he meant to make “an attempt to develop those branches 
of physics which are of fundamental importance to many of the 
operations concerned in industrial chemistry.” 

The expectation raised that here is a practical book for 
industrial chemists and chemical engineers is realized only for 
a very small part, as the author does not sufficiently bring out 
how and when the knowledge contained in the book can be 
applied in this industry. 

In 700 pages the au hor discusses a fair section of mathe- 
matics, almost the entire physical chemistry, the physics of 
properties of matter and a part of chemical engineering. In 
consequence, the book has become a collection of short essays 
which on the one hand are too elementary for allowing the 
book to be classed among the average scientific publications, 
while on the other hand they are written in such a compact 
style as to be only readable for the expert. For these reasons 
it might be doubted whether the book will be a proper guide 
for the student. 

The book is composed of four sections: 

Part “Mathematical Introduction.” 
with the entire mathematical apparatus which in the author's 
opinion is needed by the chemical engineer; standards and 
measurements, vector and tensor quantities, mechanical ma- 
thematies (i.a. nomography, graphical representation, Heavisi- 
de, operational calculus) dimensional analysis and statistical 


This part deals 


methods. 

Part II: “ Properties of matter.” This part, which more or 
less resembles a textbook of physical chemistry (however, 
without bearing comparison with other works on this matter), 
also deals with rheology. 

Part III treats of movement of matter and deals with: type 
of motion; wave motion; diffusion; diffusion of mass, moment - 
um, energy and electric charge; simultaneous diffusion of two 
or more quantities and transflux phenomena. 


Part IV entitled “ Scientific Instruments,” starts with an 
introduction on measuring, recording and controlling instru- 
ments, whereafter the author discusses the technique of me- 
asuring length, surface, volume, angles, velocity, mass, pressu- 
re, Viscosity, elasticity, fluid flow, optical and electrical pro- 
perties and temperature. This part is written in an elementary 
style but gives a clear explanation of the underlying theory. 

The entire book gives evidence of the author's very exten- 
sive knowledge of the matter. It is properly formulated. The 
general arrangement of the book is to be highly praised, espe- 
cially in Part III where heat and mass transfer, fluid flow, 
electricity etc. are treated from a uniform point of view. It 
is to be regretted that as a starting point the author did not 
select the general differential equation which underlies all 
these processes. 

Part LV deserves much appreciation, especially the chapter 
on temperature measurements. 

Remarkable features of this book are the far too succinct 
discussion of the systems of units and the little attention paid 
to electronics and modern spectrometry. 


Book reviews 


Chemical 
Engineering Science 


The book can be consulted as an encyclopaedia; when used 
in this manner it is always interesting to see what the author 
has observed about a given subject and for this purpose it 
can certainly be recommended. 

A decided handicap is that the book, which was written 
between 1942 and 1945, did not appear in print until 1950. 

It is to be hoped that in the next edition the number of 
subjects will be greatly reduced and that the matter will be 
dealt with in a manner which is more directly aimed at in- 


dustrial application. W. vaN KreveLen 


Dr. Ing. Destillier- und Rektilfiziertechnik. 
Springer-Verlag, Berlin-Gottingen-Heidelberg, 1950. Zweite 
Auflage, XVI 465 pp. mit 204 Abbildungen im Text und 
23 Kurventafeln, 16 » paper DM 45,—, bound 
DM 49.50 


24 em. Price: 


Good wine needs no bush: the publication of the second 
edition of Prof. Kirsenpaums well-known book is to be warmly 
welcomed. As compared with the first edition (and its trans- 
lation) the edition abundance of new 
material, the greater part of which has been borrowed from 


second contains an 
the American literature. 

The main division of the book is as follows: theory (units; 
vapour-liquid equilibria, partial condensation and heat of 
vaporization of mixtures); simple distillation; the rectifying 
column (basic principles, construction, caleulation of number 
of theoretical plates, minimum reflux, heat consumption) 
continuous rectification (principles, calculation, special cases, 
heat losses); the enthalpy-concentration diagram; multi- 
component separation (ideal three- and multicomponent mix- 
tures); design of plate columns (influence of special construct. 
ion, comparison between bubble-cap and perforated plates, 
influence of fluid flow); design of packed columns (influence 
of packing construction, liquid flow, concentration ete. pressu- 
re losses); liquid hold up in columns; heat exchangers and ope- 
rational control (brief). The author has aimed at presenting 
the distillation and rectifying technique as a whole by omitting 
the molecular distillation; extractive distillation is dealt with 
in a cursory manner. 

The argumentation and explanations are exact and de- 
tailed; much attention is given to the problems of heat trans. 
Of great 
value are the discussions of the rectification of multicompo- 


fer, heat balance and to joint heat and mass transfer. 


nent systems, the comparison of batch and continuous di- 
stillations and the data on packed columns. 

As regards printing and illustration the book is well got-up, 
plan and finish being typically European. Units and symbols 
will be uncommon to the Anglo-Saxon reader. 


The large number of tables and graphs in which as many as 
162 binary and 11 ternary vapour-liquid equilibria are fully 
described are most valuable. The most important physico- 
chemical constants of 114 simple substances occurring in these 
systems are given. This would be sufficient to make the book 
a desirable possession. 


Briefly summarized it may be said that this book will be 
a Valuable guide for graduate and postgraduate study while 
experts in this field will find it an important source of in- 


formation. D. W. van Kreve ten 


